vl/ 


DEVELOPMENT  AND  TESTING  OF  A  MULTIPURPOSE 
HYDROLOGIC  YIELD  MODEL 

by 

JUDITH  MARIE  HAYDEN 
B.5.,  Oklahoma  State  University,  1977 


A  MASTER'S  THESIS 


submitted  In  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  cf  Civil  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 
1979 


Ape roved  by: 


~ 


#^  y^^^ts 


Major  Professor^ 


LD 

17 

h: 
x 

ACKNOWLEDGMENTS 
The  author  wishes  to  express  her  gratitude  and  appreciation  to  the 
following  individuals:   to  her  major  advisor,  Dr.  J.  J.  Zovne,  for  his 
encouragement  and  consultation  throughout  this  research  study;  to  her 
committee  members,  Drs.  J.  K.  Koelliker  and  H.  L.  Manges,  for  their 
encouragement  and  assistance  during  this  study;  to  her  husband,  Myron, 
for  his  encouragement  and  understanding  throughout  her  graduate  studies; 
to  her  parents  for  being  dedicated  individuals  from  which  to  learn;  to 
Mrs.  Joan  Edwards  for  typing  the  manuscript;  and  to  the  Engineering 
Experiment  Station,  Agricultural  Experiment  Station,  Department  of 
Civil  Engineering,  Kansas  Water  Resources  Research  Institute  and  the 
State  of  Kansas  Water  Resources  Board  for  providing  financial  support 
during  this  research. 


TABLE  OF  CONTENTS 

Chapter  Page 

I.   INTRODUCTION  1 

Statement  of  Problem   1 

II.   REVIEW  OF  PREVIOUS  PROGRAM  PARAMETERS   3 

Evapotranspiration   3 

Surface  Runoff   7 

Moisture  Redistribution  3 

III.   DEVELOPMENT  OF  NEW  METHODOLOGICAL  CONCEPTS  11 

Weather  Modification   11 

Moisture  Redistribution  16 

Municipal  Wastewater  Effluents   19 

IV.   DESCRIPTION  OF  COMPUTER  PROGRAM   25 

Purpose 25 

Organization   26 

Capabilities   39 

Data  Input 41 

Data  Output 50 

V.   RESULTS  AND  DISCUSSION 53 

Weather  Modification   53 

Municipal  Wastewater  Effluents   62 

VI.   CONCLUSIONS  AND  RECOMMENDATIONS   66 

Conclusions 66 

Recommendations  67 

REFERENCES 68 

APPENDIX  A  -  LIST  OF  SYMBOLS 70 

APPENDIX  3  -  REFERENCE  IDENTIFIERS  73 

APPENDIX  C  -  LISTING  OF  COMPUTER  PROGRAM 78 

APPENDIX  D  -  INPUT  FOR  COMPUTER  ANALYSIS 104 

APPENDIX  E  -  OUTPUT  OF  COMPUTER  ANALYSIS  109 


ixi 


LIST  OF  TABLES 

Table  Page 

I.   Precipitation  Enhancement  Models  15 

II.   Guide  for  Data  Input 42 

III.   Summary  of  Results  for  Precipitation  and  Runoff   54 

IV.  Summary  of  Results  for  Recharge  and  Evapotranspiration  .  .  61 

V.  Summary  of  Results  for  a  Disposal  Control  Facility  ....  64 

VI.  Runoff  Curve  Numbers  for  Antecedent  Soil  Moisture 

Condition  II 74 

VII.   Soil  Classification  Used  for  Irrigation  Design 75 

VIII.   Crop  Codes  Used  in  the  Computer  Program 77 

IX.   Area  Cedes  Used  in  Che  Computer  Program 77 

X.   INPUT  DATA  FOR  NAMELIST  ALPHA 105 

XI.   INPUT  DATA  FOR  NAMELIST  BETA 106 

XII.   INPUT  DATA  FOR  NAMELIST  AREAS 108 


lv 


LIST  OF  FIGURES 

Figure  Page 

1.  Schematic  of  the  Subsurface  Profile   9 

2.  Schematic  Representation  of  Moisture  Redistribution   ....  18 

3.  Generalized  Flowchart  for  Main  Program  27 

4.  Generalized  Flowchart  for  Subroutine  CROPCO   32 

5.  Generalized  Flowchart  for  Subroutine  WTRMOD   33 

6.  Generalized  Flowchart  for  Subroutine  SNOWRT   35 

7.  Generalized  Flowchart  for  Subroutine  DARCRT   36 

8.  Generalized  Flowchart  for  Subroutine  STORAG   38 

9.  Natural  and  Enhanced  Annual  Precipitation  for 

Western  Kansas   57 

10.  Natural  and  Enhanced  Annual  Precipitation  for 

Central  Kansas   58 

11.  Natural  and  Enhanced  Annual  Precipitation  for 

Eastern  Kansas   59 


v 


CHAPTER  I 

INTRODUCTION 

A  dependable  supply  of  water  for  domestic,  agricultural  and  industrial 
use  in  the  United  States  has  increased  in  importance  during  the  past  few 
years.   Although  water  is  fundamentally  a  renewable  resource  in  some  regions, 
it  is  becoming  increasingly  scarce.   For  example,  the  water  levels  in  many 
reservoirs  are  falling  and  groundwater  tables  are  declining  to  dangerously 
low  levels  in  many  areas,  thus  resulting  in  periodic  restrictions  placed 
upon  its  use.   In  an  effort  to  remedy  these  conditions  extensive  research 
is  being  conducted  on  the  technological  possibilities  of  providing  suffi- 
cient water  to  meet  man's  increasing  demands.   These  technological  possibili- 
ties include  providing  new  supplies  of  fresh  water  while  increasing  the 
efficiency  of  available  supplies. 

Statement  of  Problem 

A  possible  method  for  increasing  the  water  supply  within  a  particular 
locality  is  by  weather  modification.   Weather  modification  is  an  attempt 
to  alter  the  atmospheric  conditions,  thereby  providing  a  more  favorable 
distribution  of  precipitation.   During  the  past  thirty  years  the  main 
emphasis  of  weather  modification  has  been  to  increase  localized  precipita- 
tion.  Although  considerable  research  has  been  conducted  in  the  past  two 
decades  on  cloud  seeding  technology,  the  effect  on  the  hydrological 
environment  remains  uncertain. 


In  addition  to  methods  of  providing  new  supplies  of  water  by  weather 
modification,  methods  to  increase  the  efficiency  of  available  water  supplies 
have  also  received  considerable  attention.   The  method  of  optimizing 
available  water  supplies  considered  in  this  study  is  the  utilization  of 
effluent  from  wastewater-treatment  plants  for  irrigation  purposes.   This 
method  involves  the  controlled  discharge  of  wastewater  to  the  land  for  the 
purpose  of  promoting  plant  growth.   A  major  factor  which  must  be  considered 
in  land  application  of  wastewater  effluent  is  the  design  of  a  storage 
facility  for  wastewater  control.   The  design  of  such  a  facility  requires 
the  evaluation  of  the  relative  effect  of  various  conditions  on  the 
minimum  containment  volume  necessary  to  store  and  control  the  wastewater 
effluent. 

This  research  study  consisted  of  the  development  of  a  new  methodo- 
logical approach  to  the  previously  described  problems.   The  new 
methodology  was  incorporated  into  a  computer  program  initially  developed 
by  Bean  (1) .   The  final  computer  program  considers  selective  land 
disposal  and  evaporative  techniques  to  enhance  agricultural  production 
under  various  geographic  and  climatic  conditions.   The  computer  program 
can  be  used  to  predict  the  optimized  design  for  a  disposal  control 
facility  in  addition  to  evaluating  the  hydrological  effects  resulting  from 
the  implementation  of  a  rainfall  alteration  scheme. 


CHAPTER  II 

REVIEW  OF  PREVIOUS  PROGRAM  PARAMETERS 

The  purpose  of  this  Chapter  is  to  review  the  basic  theories  and  con- 
cepts utilized  in  the  computer  program  developed  previous  to  this  study 
which  established  feedlot  runoff  control  and  disposal  guidelines  (1) . 
Presented  in  the  first  section  are  the  concepts  and  principles  of  evapo- 
transpiration.   The  second  section  reviews  the  hydrological  aspects  of 
infiltration,  interception  and  surface  runoff.   Presented  in  the  final 
section  are  the  concepts  of  moisture  redistribution. 

Evapo transpiration 

The  concept  of  evapotranspiration  was  incorporated  into  the  program  to 
evaluate  the  change  in  the  soil-moisture  regime  as  a  result  of  a  change  in 
the  climatic  environment.   Evapotranspiration  as  used  in  this  study  is 
defined  as  the  rate  of  moisture  loss  to  the  atmosphere  as  a  result  of  the 
evaporation  of  water  from  plants  and/or  a  specified  surface  under  ambient 
climatological  conditions.   It  should  be  noted  that  a  difference  exists 
between  potential  and  actual  evapotranspiration  for  any  given  site.   The 
potential  evapotranspiration,  more  correctly  termed  the  evapotranspiration 
potential,  is  the  maximum  rate  of  moisture  lost  under  ideal  conditions. 
Therefore,  the  evapotranspiration  potential  is  the  upper  boundary  limit  on 
the  rate  of  moisture  loss.   This  fact  necessitates  that  the  actual  value 
of  the  evapotranspiration  can  never  be  greater  than  the  evapotranspiration 


potential  and  is  usually  appreciably  lower.   The  evapotranspiration  potential 
as  well  as  the  actual  evapotranspiration  for  any  crop  and  soil  condition 
can  only  be  estimated  by  indirect  means. 

A  common  method  for  computing  evapotranspiration  is  the  Penman 
combination  equation.   This  equation  is  capable  of  considering  the  effects 
of  various  climatic  conditions  on  the  rate  of  evaporation  likely  to  occur 
from  lakes,  bare  soil,  and  vegetative  surfaces.   The  Penman  equation 
incorporates  both  energy  balance  and  aerodynamic  transport  theories  into 
a  single  mathematical  expression.   The  mathematical  model,  based  on  the 
Penman  equation  as  developed  by  Bean  (1) ,  for  determining  the  daily  evapo- 
transpiration is  shown  in  Equation  (1) . 

ETP  =  0.039  Ta°'673  [ (1-r)  Ra  (0.22  +  0.54  PSUNS)  -  (   . 

2.010  x  10"9  T4  (0.98  -  a-b  /ES  x  RHD)x(0.1  +  0.9  PSUNS)] 
+  (1  -  0.039  Ta°'673)  x  0.26(e  +  0.01  WD)  (ES  -  ES  x  RED) 


where 


ETP  =  evapotranspiration  potential,  in  inches 

Ta  =  mean  daily  air  temperature,  in  F 

T  =  mean  daily  air  temperature,  in  K 

Ra  =  solar  radiation,  in  mm  of  water 

PSUNS  =  percent  sunshine,  in  %/100 

ES  =  saturated  vapor  pressure  at  Ta,  in  millibars 

RHD  =  relative  humidity,  in  %/i00 

WD  =  wind  speed,  in  miles/day 

a,b  -  geographic  location  constants 

r  =  reflectance  coefficient  (albedo) 

e  =  mass  transfer  coefficient 


5 
The  variables  Ra,  PSUNS ,  RHD  and  WVD  for  a  specific  geographic 
location  are  extrapolated  monthly  averages  from  first-order  weather  sta- 
tion data  while  T  and  Ta  are  computed  from  daily  climatological  records. 
The  value  of  the  reflectance  coefficient  is  dependent  on  the  type  of 
surface  cover.   The  evaporation  potential  is  equivalent  to  evaporation 
from  a  free  water  surface  when  r  is  equal  to  0.05.   For  green  crops  r 
varies  from  0.20  to  0.25.   The  geographic  location  constants  can  be 
obtained  by  the  methods  described  by  Zovne  and  Koeiliker  (2) .   The  saturated 
vapor  pressure,  according  to  Linsley  (3),  can  be  computed  by  the  use  of 
Equation  (2) . 

ES  -  33.9  [(0.00738  +  0.8072)8  -  0.000019  |  l.STa  +  48  | 

(2) 
+  0.001316] 

where 

ES  and  Ta  have  been  previously  defined. 

The  actual  rate  of  evapotranspiration  for  either  bare  or  vegetated 
soil  conditions  is  affected  by  the  type  of  soil  and  crop  present.   Evapo- 
transpiration is  believed  to  occur  in  two  stages  (4).   The  first  stage, 
known  as  che  constant  rate  stage,  occurs  when  the  available  soil  moisture 
is  greater  than  thirty  percent  of  the  maximum  amount  of  available  soil 
moisture.   Hillel  (4)  defines  available  soil  moisture  as  the  difference 
between  the  in  situ  soil  moisture  concent  and  the  permanent  wilting  point. 
The  maxirium  available  soil  moist  .ire  is  defined  as  the  difference  between 
the  field  capacity  and  the  permanent  wilting  point.   To  calculate  the  first 
stage  evaporation,  Equation  (1)  is  used  for  a  bare  soil  condition  with  r 
equal  to  0.20. 
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When  the  available  soil  moisture  is  less  than  thirty  percent  of  the 
maximum  available  soil  moisture,  second  stage  evapotranspiration  occurs. 
To  calculate  the  evaporation  from  bare  soil  for  this  stage  Equation  (3) 
can  be  used. 

E  -  c'  t2   -  c*  (t  -  l)3*  (3) 

s 

where 

E  =  soil  evaporation,  in  inches 

c'  =  hydraulic  coefficient  of  the  soil,  in  inches/day 
t  =  time  after  the  first  stage  evaporation,  in  days 
If  a  vegetated  surface  exists,  then  Equation  (4)  should  be  used  to  cal- 
culate the  evapotranspiration. 

AET  =  ETP  x  k  x  9   /  .3  9  (4) 

a       max 

where 

AET  =  actual  evapotranspiration,  in  inches 

k  =  crop  consumptive  use  coefficient 

9  =  available  soil  moisture,  in  inches 
a 

9    =  maximum  available  soil  moisture,  in  inches 
max 

ETP  =  same  as  defined  previously 
When  the  soil  moisture  content  is  at  the  permanent  wilting  point,  evapo- 
transpiration does  not  occur.   Equation  (4)  therefore  assumes  a  linearly 

decreasing  AET  from  ETP  to  zero  from  0.3  9    to  wilting  point.   The  crop 
b  max 

consumptive  use  coefficient  (K)  incorporates  the  Blaney-Criddle  Method  (5) 
for  evaluating  actual  evapotranspiration.   A  daily  k  value  is  computed  as 
described  in  (2) . 


Surface  Runoff 

A  method  for  estimating  surface  runoff  which  has  received  wide 
acceptance  has  been  developed  by  the  Soil  Conservation  Service  (SCS)  (6) . 
This  method  assumes  that  there  is  an  initial  abstraction  (IA)  of  rainfall 
prior  to  the  occurrence  of  any  surface  runoff.   The  initial  abstraction 
consists  of  infiltration,  interception  and  surface  storage.   Incorporated 
into  the  method  is  the  assumption  that  IA  is  equal  to  0.2  of  the  maximum 
potential  dif ference(s)  between  precipitation  and  runoff.   Equation  (5) 
can  be  used  to  calculate  this  maximum  difference. 

S  -  32pfi  -  10  (5) 

where 

N  =  runoff  curve  number 

S  =  maximum  potential  difference,  in  inches 
As  can  be  seen  in  this  Equation,  3  is  dependent  upon  the  antecedent  moisture 
conditions,  soil  conditions,  land  use,  and  conservation  practices  through 
the  use  of  a  value  (N)  obtained  from  generalized  runoff  curves.   Values 
for  N  are  listed  in  Table VI  located  in  Appendix  B. 

Upon  evaluating  the  maximum  potential  difference,  the  magnitude  of 

surface  runoff  can  be  estimated  by  the  use  of  Equation  (6) . 

2 
) 

P  + 
where 

Q  =  surface  runoff,  in  inches 
P  =  precipitation,  in  inches 
S  =  as  defined  previously 


.    (P  -  0.2  S)  ,,, 

Q  =  5  I  n  q  c —  (6) 
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This  equation  is  applicable  only  when  the  amount  of  precipitation  is 
greater  than  twenty  percent  of  S.   When  precipitation  is  less  than  twenty 
percent  of  S,  surface  runoff  is  assumed  not  to  occur. 

To  evaluate  the  magnitude  of  the  initial  abstraction  believed  lost 
as  a  result  of  interception  Bean  (1)  suggested  using  an  interception- 
storage  capacity  of  0.1  inch  which  is  evaporated  at  the  free  potential 
rate.   By  this  method  the  interception-storage  requirement  must  be  satisfied 
prior  to  the  evaluation  of  any  infiltration  or  runoff.   The  remainder  of 
the  potential  difference  which  is  not  lost  to  evaporation  or  interception 
is  considered  as  infiltration. 

Moisture  Redistribution 

This  section  of  the  Chapter  discusses  the  concepts  used  by  Bean  (1) 
for  the  redistribution  of  moisture  within  the  soil  profile.   The  soil 
profile  was  divided  into  two  layers  with  the  first  layer  being  one  foot 
in  thickness  and  the  second  layer  being  three  feet  thick.   These  two 
layers  were  believed  to  adequately  describe  the  potential  root  zone 
within  the  subsurface  soil  profile.   This  assumption  is  presented  con- 
ceptually in  Figure  1. 

Based  on  this  assumption  infiltration  which  resulted  from  either 
rainfall  or  irrigation  was  assumed  to  fill  each  layer  to  ninety  percent 
of  the  saturation  value  with  any  excess  moisture  passing  below  four  feet 
to  be  used  as  groundwater  recharge.   When  the  time  between  any  successive 
events  which  could  result  in  infiltration  exceeded  two  days,  the  soil 
moisture  of  the  upper  zone  was  decreased  to  field  capacity  with  the  excess 
moisture  redistributed  to  the  lower  zone.   The  lower  zone  was  then  decreased 
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Figure  1.   Schematic  of  the  Subsurface  Profile 
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to  ninety  percent  field  capacity.   The  remaining  moisture  was  considered 
to  percolate  below  the  two  layers  as  a  source  for  groundwater  recharge. 
This  method  for  the  redistribution  of  moisture  has  produced  reasonable 
estimates  of  the  vertical  movement  of  infiltrating  water,  although  it 
neglects  upward  movement  of  sub-soil  water. 


CHAPTER  III 

DEVELOPMENT  OF  NEW  METHODOLOGICAL  CONCEPTS 

Presented  in  this  Chapter  are  the  basic  theories  utilized  in  the  develop- 
ment of  new  methodological  concepts  incorporated  in  the  computer  program  (1) . 
The  Chapter  is  divided  into  three  sections.   Tha  first  section  reviews  the 
predictive  techniques  to  evaluate  the  effect  of  weather  modification  on 
precipitation.   The  concepts  and  principles  of  moisture  redistribution  within 
the  upper  subsurface  stratum  are  presented  next.   Presented  in  the  final 
section  are  the  current  theories  and  research  which  relate  to  the  use  cf 
municipal  waste  as  a  controllable  source  of  water  for  irrigation. 

Weather  Modification 

Weather  modification  encompasses  all  aspects  of  weather  and  climate 
changes  brought  about  by  the  activities  of  man.   This  section  of  the 
Chapter  presents  the  current  theories  of  weather  modification  which  could 
be  used  to  increase  the  amount  of  precipitation  within  a  particular  locality. 
The  principal  technique  used  is  cloud  seeding  from  airborne  or  ground- 
based  generators. 

Two  methods  whereby  precipitation  could  be  increased  artificially 
by  cloud  seeding  have  been  suggested  by  Hobbs  (7) .   The  first  method  is 
the  introduction  of  large  hygroscopic  particles  or  liquid  water  droplets 
into  a  warm  cloud  formation  which  lies  completely  below  the  level  of  the 
0°C  isotherm.   The  artificial  nuclei  initiate  the  coalescence  mechanism 
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of  rainfall  production  by  falling  through  the  cloud  which  results  in  the 
collision  and  coalescence  of  the  smaller  droplets  lying  in  their  paths. 
In  this  way  the  larger  droplets  tend  to  increase  in  size  rapidly.  Upon 
emerging  from  the  base  of  the  cloud  the  larger  droplets  should  be  large 
enough  to  survive  evaporation  in  the  dry  air  and  reach  the  ground  as 
raindrops. 

The  second  method  to  increase  precipitation  is  to  introduce  artificial 
ice  nuclei  in  a  cold  cloud.   Precipitation  from  a  cold  cloud  is  dependent 
upon  the  coexistence  of  supercooled  water  droplets  and  ice  particles  in  the 
regions  of  a  cloud  above  the  0°C  isotherm.   Since  the  concentration  of 
natural  ice  nuclei  in  the  atmosphere  is  small,  cloud  seeding  with  artifi- 
cial ice  nuclei  could  result  in  the  rapid  growth  of  the  ice  crystals  by 
condensation  of  the  ambient  moisture.   The  ice  crystals  then  combine  to 
form  snowflakes  which  fall  to  the  0  C  isotherm,  and  subsequently  melt  to 
form  raindrops  as  they  reach  the  ground. 

During  the  past  thirty  years  experimental  research  has  been  conducted 
concerning  the  effects  of  cloud  seeding  on  the  amount  of  precipitation. 
The  first  extensive  experiment  conducted  on  the  seeding  of  warm  clouds  was 
carried  out  in  1952.   Warm  clouds  were  seeded  with  small  particles  of 
sodium  chloride  in  East  Africa  by  exploding  bombs  impregnated  with  sodium 
chloride  just  above  the  base  of  the  cloud  (7) .   On  some  occasions  the 
clouds  in  which  the  sodium  chloride  bombs  were  exploded  were  noted  to 
disperse  rapidly.   However,  during  the  seeding  period  of  one  month  the 
total  rainfall  in  an  area  10  miles  downwind  had  increased  six  inches  over 
a  similar  period  in  this  area  of  nonseeded  days. 
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During  this  same  period  of  time  Bowen  (7)  experimented  in  Australia 
by  seeding  warm  clouds  directly  with  water  droplets.   In  this  experiment, 
the  behavior  of  the  seeded  clouds  were  compared  with  that  of  similar 
nonseeded  formations  and  an  increase  in  precipitation  was  observed.   More 
recently,  the  University  of  Chicago  conducted  research  by  seeding  warm 
cumulus  clouds  in  the  Caribbean  with  water  droplets.   An  increase  in 
clouds  developing  precipitation  was  noted. 

The  first  experimental  research  conducted  on  seeding  cold  clouds 
occurred  in  1946  by  Vincent  Schaef er  and  Irving  Langmuir  (8) .   Schaef er 
and  Langmuir  observed  that  dry  ice,  when  dropped  into  a  cloud  of  super- 
cooled droplets  under  laboratory  conditions,  produced  numerous  ice  crystals. 
Consequent  field  testing  verified  dry  ice  to  be  an  effective  means  to 
nucleate  large  areas  of  supercooled  clouds. 

Following  this  demonstration,  Bernard  Vonnegut  (8) ,  who  was  also 
working  with  Schaef er  and  Langmuir,  began  to  search  for  other  substances 
having  the  same  crystal  properties  as  ice.   Vonnegut  found  that  a  silver 
iodide  solution  was  an  effective  ice  nucleator  and  could  be  applied  from  a 
ground-based  generator,  thereby  reducing  the  cost  of  modifying  cold  clouds. 
Thus,  a  silver  iodide  solution  became  the  source  for  consequent  cloud 
seeding  research. 

In  1966  the  National  Academy  of  Sciences  Panel  on  Weather  and  Climate 
Modification  analyzed  the  results  of  research  projects  involving  cloud 
seeding  with  ground-based  silver  iodide  generators.   The  panel  concluded 
that  the  artificial  seeding  of  clouds  apparently  could  modify  cloud 
structure  and  increase  precipitation  by  an  average  of  approximately  ten 
percent.   However,  it  remained  uncertain  whether  it  is  possible  to  increase 
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precipitation  over  large  areas  for  extended  periods  of  time  by  this 
weather  modification  technique.   Many  more  experiments  need  to  be  con- 
ducted before  definitive  conclusions  can  be  drawn. 

To  estimate  the  potential  effects  of  weather  modification  on  the 
environment  it  was  assumed  in  this  study  that  cloud  seeding  can  produce 
changes  in  precipitation  and  that  this  change  will  occur  in  a  specified 
manner.   Based  on  the  research  conducted  by  Changnon  and  Huff  (9)  the 
effects  of  a  seeding  program  will  produce  a  seventy-five  percent  increase 
in  precipitation  of  normally  small  precipitaticnal  clouds.   They  further 
indicate  that  clouds  producing  normally  large  amounts  of  precipitation 
do  not  benefit  from  cloud  seeding  and,  in  some  cases,  seeding  can  reduce 
the  natural  efficiency  of  a  cloud  by  ten  percent.   Factors  hindering  the 
effectiveness  of  a  cloud  seeding  program  are  incorporated  into  the 
precipitation  alteration  scheme  developed  during  this  research  study. 

Model  3-ZERO  which  is  presented  in  Table  1  is  based  on  the  results 
of  Changnon  and  Huff's  research  (10).   The  model  alters  the  daily 
rainfalls  in  Kansas  by  the  percentages  indicated  during  the  months 
March  through  September.   In  the  model  a  discontinuity  occurs  at  the 
24-hour  precipitation  value  of  one  inch.   To  correct  this  anomoly  and 
also  to  provide  a  precipitation  enhancement  range,  Model  B-ZERO  MODIFIED 
was  developed  (11).   Model  B-ZERO  MODIFIED,  also  presented  in  Table  I, 
assumes  that  there  is  no  increase  nor  decrease  to  daily  rainfalls  of 
greater  than  one  inch.   The  effects  of  these  two  precipitation  alteration 
models  on  the  hydrological  aspects  in  Kansas  are  discussed  in  detail  in 
Chapter  V. 
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Moisture  Redistribution 

The  process  of  moisture  redistribution  within  a  soil  profile  primarily 
occurs  under  unsaturated  conditions.   The  unsaturated  flow  process  is 
difficult  to  describe  quantitatively  because  of  the  complex  relationships 
between  water  content,  hydraulic  conductivity,  and  hydraulic  gradient. 
Hillel  (4)  suggested  that  an  accepted  method  for  estimating  the  movement 
of  water  one-dimensionally  is  by  use  of  the  Darcy  equation,  represented 
mathematically  by  Equation  (7) . 

q'=  K  At  i  (7) 

where 

q'  =  volume  of  water  movement  per  unit  area,  in  inches 

K  =  hydraulic  conductivity,  in  inches/day 

At  =  time  increment,  in  days 

i  =  hydraulic  gradient,  in  feet /feet 

The  hydraulic  conductivity  in  Equation  (7)  is  a  function  of  the 
water  content  of  the  soil.   The  value  of  the  hydraulic  conductivity  is  at 
its  greatest  when  the  soil  is  saturated,  and  decreases  steeply  with 
decreasing  water  content.   Since  the  value  of  soil  moisture  is  a  continuous 
variable  over  time  within  the  computer  program,  unsaturated  hydraulic 
conductivity-soil  moisture  content  relationships  for  the  in  situ  soil  are 
required  to  apply  Equation  (7) .   These  relationships  can  be  obtained  from 
laboratory  experimentation. 

The  hydraulic  gradient  in  Equation  (7)  is  the  sum  of  the  gravitational 
gradient,  directly  dependent  on  elevation,  and  the  soil  water  potential. 
This  is  expressed  mathematically  by  Equation  (8) . 
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Az  +  Ah 

1  =  : (8) 

Az 

where 

Az  =  change  in  elevation,  in  feet 

Ah  =  change  in  soil  water  potential  over  distance  Az,  in  feet 

i  is  the  same  as  previously  defined 
The  soil  water  potential  is  assumed  to  be  a  positive  value  when  the  water 
movement  is  in  the  downward  direction  indicating  a  decrease  in  soil  water 
potential  with  a  decrease  in  elevation.   A  negative  soil  water  potential, 
also  referred  to  in  literature  as  a  capillary  potential  or  a  soil  water 
suction,  occurs  when  the  decrease  in  soil  water  potential  is  in  the  upward 
direction  causing  an  upward  movement  of  flow.   To  apply  the  value  of  the 
soil  water  potential  within  the  computer  program  to  Equation  (7)  laboratory 
experimentation  relating  soil  moisture  and  the  soil  water  potential  is 
required. 

To  evaluate  the  redistribution  of  moisture  using  Darcy's  equation 
for  one-dimensional  flow,  che  subsurface  profile  illustrated  in  Figure  2 
was  incorporated  conceptually  into  the  computer  program.   As  illustrated, 
the  subsurface  profile  was  divided  into  three  layers  with  the  upper  layer 
being  one  foot  in  thickness,  the  middle  layer  being  three  feet  in  thick- 
ness, and  the  lower  layer  being  two  feet  thick.   The  three  layers  were 
believed  to  adequately  describe  the  movement  of  moisture  within  the  poten- 
tial root  zone  of  the  subsurface  profile. 

An  infiltrating  event  resulting  from  either  rainfall  or  irrigation  is 
distributed  within  the  soil  profile  by  filling  each  successive  layer  to 
field  capacity.   All  excess  moisture  passing  below  six  feet  is  assumed  to 
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Figure  2.   Schematic  Representation  of  Moisture  Redistribution 
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be  eventually  used  as  groundwater  recharge.   The  moisture  within  each 
layer  is  then  redistributed  by  applying  the  Darcy  equation  for  one- 
dimensional  flow.   The  time  increment  it,  suggested  by  Neibling,  et  al. 
(12),  is  assumed  to  be  equal  to  one  day  when  no  water  infiltrates  the 
soil  surface  on  any  given  day  and  equal  to  1/6  day  when  an  infiltration 
event  occurs. 

The  use  of  the  Darcy  equation  for  the  one- dimensional  unsaturated 
redistribution  of  moisture  requires  the  use  of  soil  water  potential- 
soil  moisture  and  unsaturated  hydraulic  conductivity-soil  moisture 
relationships  obtainable  from  laboratory  testing  on  undisturbed  soil 
samples,  as  described  by  Neibling,  et  al.  (12). 

Municipal  Wastewater  Effluents 

The  application  of  municipal  wastewater  effluents  on  agricultural 
land  has  received  considerable  attention  during  recent  years.   With  the 
enactment  of  the  Federal  Water  Pollution  Control  Amendments  cf  1972 
requiring  that  the  discharge  of  pollutants  into  waterways  cease  by  1935, 
the  irrigation  method  of  land  application  has  been  recognized  as  an 
alternate  means  to  treat  wastewater.   The  irrigation  process  is  used 
primarily  to  maximize  crop  production  while  the  disposing  of  municipal 
wastewater  is  a  secondary  benefit. 

There  are  many  factors  affecting  the  use  of  wastewater  effluent 
as  a  source  of  water  for  irrigation.   The  first  factor  to  be  considered 
is  the  minimum  level  of  wastewater  treatment  required  prior  to  land 
application.   Wastewater  effluent  from  primary  and/or  secondary 
treatment  plants  has  been  used  for  irrigation  in  many  states.   As 
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an  example,   California  allows  the  use  of  primary  effluent  for 
surface  irrigation  of  orchards,  vineyards,  and  fodder,  fiber  and  seed 
crops  (13).   Generally,  the  public  health  agencies  of  each  state  place 
limitations  on  the  quality  of  municipal  wastewater  that  can  be  used  for 
irrigation. 

The  next  factor  to  be  considered  is  site  selection.   Hall,  Wilding 
and  Erickson  (14)  suggest  that  the  site  selection  criteria  be  based  upon 
three  interrelated  parameters.   These  parameters  include  general  topographic 
features,  geologic  characteristics,  and  in  situ  properties  of  the  soil. 
An  evaluation  of  the  proposed  site  includes  both  a  careful  inspection  of 
the  site  and  preliminary  tests  conducted  on  the  soils. 

The  ideal  topography  for  utilization  of  wastewater  by  irrigation  is 
directly  dependent  upon  the  type  of  distribution  system  to  be  used.   For 
example,  the  irrigation  of  cultivated  crop  lands  using  the  sprinkler 
technique  restricts  the  slopes  to  fifteen  percent  or  less.   However,  slopes 
of  up  to  thirty  percent  are  permissible  if  the  sprinkler  method  of  irriga- 
tion is  to  be  applied  where  noncultivated  crops  are  grown.   Surface 
irrigation  requires  relatively  flat  land  although  slopes  as  steep  as  five 
percent  could  be  used  with  contour  furrows.   Increasing  slopes  generally 
result  in  a  reduction  of  the  infiltration  rates  of  most  soils  thus  causing 
excessive  runoff  and  soil  erosion.   The  cost  effectiveness  of  various 
irrigation  practices  is  a  major  constraint  in  selecting  the  topography 
best  suited  for  application  of  the  wastewater  effluent. 

The  suitability  of  a  site  is  also  a  function  of  the  physical, 
chemical  and  mineralogical  properties  of  the  soil.   These  properties 
influence  the  infiltration  rate,  moisture  holding  capacity  and  any 
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absorption  reactions  with  the  waste  components.   In  land  application  of 
municipal  wastewater,  the  soil  functions  both  as  a  natural  filter  and 
as  a  medium  for  any  biological  and/or  chemical  reactions,  thus  providing 
for  the  natural  treatment  of  wastewater.   An  ideal  situation  would  be  to 
have  a  soil  with  the  ability  to  assimilate  large  quantities  of  dissolved 
solids  while  remaining  highly  permeable,  however  such  a  soil  does  not 
exist.   Therefore,  the  selection  of  the  most  advantageous  soil  is 
dependent  upon  other  constraints. 

The  storage  requirement  for  the  wastewater  effluent  is  the  third 
factor  to  be  considered.   Conventional  treatment  plants  collect,  treat, 
disinfect  and  discharge  effluent  wastewater  into  available  surface  water 
on  a  continuous  basis.   The  design  of  a  land  disposal  treatment  facility, 
like  the  conventional  plant,  must  be  operated  continuously  with  the 
exception  that  the  ability  to  discharge  effluent  is  seasonal.   The  design 
of  a  land  disposal  facility  must  include  a  storage  reservoir  capable  of 
containing  the  effluent  wastewater  for  a  period  ranging  from  three  months 
in  moderate  climates  to  seven  months  in  cold  northern  states.   This 
storage  requirement  is  directly  related  to  several  climatic  factors. 
The  primary  climatic  factors  are  temperature,  precipitation  and  snow 
cover.   The  application  of  effluent  to  the  land  is  generally  suspended 
when  the  temperature  falls  below  25  degrees  Fahrenheit  (13) .   The  maximum 
daily  precipitation  allowed  before  the  application  of  wastewater  must  be 
suspended  is  dependent  on  the  infiltration  rate  of  the  soil.   Snow  cover 
should  be  minimal  during  the  land  application  process. 

The  treatment  process  of  wastewater  by  land  application  requires  an 
extensive  monitoring  program.   The  program  should  be  designed  to  consider 
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influent,  pretreatment,  storage,  disinfection,  drainage,  seepage,  ground- 
water, surface  water,  soil  and  crop  characteristics.   Data  obtained  from 
these  monitoring  programs  is  essential  to  the  management  decisions  con- 
cerning land  application. 

In  recent  years  a  number  of  experimental  and  operational  land  appli- 
cation systems  have  been  designed  to  utilize  wastewater.   The  Michigan  State 
University  Water  Quality  Management  Project  (WQMP)  is  an  example  of  a 
research  and  development  project  which  studies  the  applications  of  waste- 
water to  agricultural  land  (15) .   This  facility,  which  was  completed  in 
1974,  utilizes  effluent  from  the  East  Lansing  sewage  treatment  plant  after 
it  has  undergone  primary  and  secondary  treatment  processes.   The  site  is 
carefully  monitored  and  controlled  to  determine  possible  changes  in 
groundwater,  soil,  or  other  variables. 

The  WQMP  facility  consists  of  four  artificial  lakes  having  a  total 
surface  area  of  forty  acres  with  an  average  depth  of  eight  feet.   The 
municipal  effluent  undergoes  chemical,  biological  and  physical  treatment 
while  it  passes  sequentially  through  the  four  lakes  over  a  period  of  30 
to  60  days.   Up  to  two  million  gallons  per  day  of  effluent  can  be  handled 
by  this  facility. 

The  facility  at  Michigan  State  University  also  includes  three  hundred 
and  twenty  acres  for  application  of  treated  wastewater  effluent  to 
agricultural  lands.   One  hundred  and  fifty  of  these  acres  have  been 
provided  with  irrigation  spray  equipment  capable  of  applying  effluent 
at  a  rate  of  two  inches  per  week  between  the  months  March  and  November. 
Both  annual  and  perennial  forage  crops  have  been  irrigated  using  this 
technique.   Although  the  project  is  still  in  the  experimental  stage,  present 
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results  indicate  the  economic  feasibility  of  applying  wastewater  effluents 
on  agricultural  land  for  the  purpose  of  enhancing  food  and  fiber  production 
for  livestock. 

The  largest  wastewater  land  treatment  facility  using  sprinkler 
irrigation  equipment  in  the  United  States  at  the  present,  time  is  the 
Muskegon  County  Wastewater  Management  System  (16) .   This  facility  was 
designed  to  transmit  municipal  and  industrial  effluents  on  a  ccunty-wide 
basis  into  a  central  treatment  plant.   At  the  central  plant  the  wastewater 
is  treated  to  comply  with  the  local,  state  and  federal  discharge  require- 
ments prior  to  land  application.   All  facets  of  the  operation  are  con- 
tinuously monitored  to  guard  against  possible  pollution  of  the  soil, 
groundwater  and  surface  waters. 

The  Muskegon  County  facility  was  designed  for  an  influent  flow  rate 
of  forty-two  million  gallons  per  day.   The  current  average  flow  is  approxi- 
mately twenty-eight  million  gallons  per  day  with  forty  percent  domestic 
wastewater  and  sixty  percent  industrial.   The  wastewater  is  initially 
treated  biologically  in  three  aerated  lagoons.   The  wastewater  is  then 
discharged  into  two  large  storage  lagoons  which  have  a  maximum  design 
storage  period  of  four  months.   The  water  quality  of  the  storage  facilities 
is  monitored  daily.   When  the  water  quality  meets  the  state  health  depart- 
ment standards,  the  treated  wastewater  is  released  from  the  storage  lagoons 
into  an  outlet  lagoon  where  it  is  disinfected  by  chlorination  prior  to 
land  application.   The  disinfected  effluent  is  then  sprayed  on  6,000  acres 
of  farmland  at  an  application  rate  of  approximately  3.5  inches  per  week. 

The  Muskegon  County  Wastewater  Management  System  has  proven  to  be  an 
efficient  means  of  treating  and  productively  disposing  of  the  municipal 
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wastewater  effluent  while  maintaining  strict  environmental  protection. 
The  elimination  of  direct  effluent  discharges  into  the  waterways  has  also 
resulted  in  improvements  in  the  water  quality  of  the  surface  water  while 
providing  conservation  of  the  groundwater. 


CHAPTER  IV 
DESCRIPTION  OF  THE  COMPUTER  PROGRAM 
Purpose 

The  purpose  of  the  computer  program  is  to  provide  a  means  of  designing 
various  disposal  control  facilities.   The  program  also  serves  as  a  means  of 
evaluating  both  selective  land  disposal  and  evaporative  techniques  to 
enhance  agricultural  production  within  Kansas.   The  general  format  of  the 
model  makes  it  applicable  over  a  wide  range  of  disposal  conditions  in 
addition  to  being  readily  adaptable  to  various  geographic  and  climatic 
regions.   The  computer  program  also  serves  as  an  effective  means  to 
evaluate  the  variation  in  hydrological  effects  which  could  result  from 
the  implementation  of  methods  to  modify  existing  weather  patterns. 

The  size  of  the  disposal  control  facility  required  is  directly 
dependent  upon  the  magnitude  of  inflow.   There  are  four  major  sources 
of  inflow  which  can  be  considered  by  this  model.   These  sources  are: 
(1)  precipitation  directly  into  the  pond,  (2)  runoff  from  the  natural 
watershed,  (3)  feedlot  runoff,  and  (4)  municipal  wastes.   Any  combination 
of  these  sources  can  be  considered  by  the  model  as  inflow  to  the  system. 
However,  outflow  from  the  system  can  occur  only  as  evaporation,  overflow 
or  induced  removal.   The  induced  removal  of  outflow  generally  consists 
of  incorporating  the  wastes  into  a  supplemental  irrigation  system.   Although 
pond  overflow  is  a  source  of  outflow  it  should  not  be  directly  considered 
in  the  design.   The  model  depicted  by  the  computer  program  simulates  field 
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conditions  by  extrapolating  the  data  obtained  from  previously  recorded 
precipitation  and  temperature  patterns.   The  model  also  considers 
variation  in  soil  properties,  meteorological  conditions,  tillage  and 
farming  practices  which  may  exist  in  the  area  under  consideration. 

Organization 

This  section  of  the  Chapter  describes  the  organizational  make-up 
of  the  computer  program.   Included  with  a  description  of  the  main  program 
are  brief  descriptions  of  each  subroutine  and  generalized  flowcharts 
illustrating  their  operational  procedures.   The  computer  program,  titled 
Kansas  State  University  Hydrologic  Yield  Model  (KSUHYDRO) ,  was  written 
in  Fortran  IV  language  and  calibrated  on  the  ITEL  AS/5  computer  located  at 
Kansas  State  University. 

KSUHYDRO  consists  of  a  main  computer  program  which  utilizes  several 
subroutines  as  a  means  of  incorporating  various  options.   The  use  of 
several  subroutines  rather  than  a  continuous  program  provides  for  a 
stable  framework  upon  which  to  base  future  modifications.   In  addition, 
this  framework  allows  for  building,  testing  and  modifying  the  program 
without  necessitating  long  periods  of  down  time. 

The  main  program,  illustrated  in  Figure  3,  controls  and  regulates 
the  various  program  components.   As  can  be  seen  in  Figure  3,  the  first 
function  of  the  main  program  is  to  read- in  the  various  information 
required.   This  information  is  read- in  on  both  data  cards  and  data 
statements.   After  the  data  has  been  entered,  properly  stored  and  cate- 
gorized the  main  program  initiates  the  sequence  of  computational  operations 
specified  by  the  model. 
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Figure  3.   Generalized  Flowchart  for  Main  Program 
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3.1 

The  computational  operations  begin  by  sizing  the  storage  facility. 
The  required  storage  facility  or  pond  is  developed  by  using  the  general 
configuration  of  a  prismatoid.   The  maximum  depth,  base  length,  base 
width,  and  side  slopes  of  the  prismatoid  are  varied  to  maximize  the 
volume  and  surface  area.   This  maximization  technique  provides  for  the 
optimization  of  required  resources. 

The  format  of  the  main  program  directs  the  computer  to  call  subroutine 
CROPCO  next.   Subroutine  CROPCO,  developed  by  Zovne,  et  al.  (2),   is  used 
to  calculate  the  various  crop  coefficients  necessary  in  the  analysis  of  the 
evapotranspiration  of  the  disposal  areas  under  consideration.   Figure  4 
shows  the  generalized  flowchart  for  this  subroutine.   Upon  returning  to 
the  main  program,  the  annual  crop  rotation  pattern  for  the  disposal  areas 
is  established. 

WTRMOD,  the  next  subroutine  called  by  the  main  program,  incorporates 
the  implementation  of  a  weather  modification  technique.   As  presented  in 
Chapter  III,  the  only  weather  modification  technique  considered  by  KSUHYDRC 
is  cloud  seeding.   The  generalized  flowchart  for  this  subroutine  is  shown 
schematically  in  Figure  5.   After  all  the  adjustments  are  made  to  the 
values  previously  recorded  for  the  daily  precipitation,  the  operational 
coding  directs  the  computer  to  return  to  the  main  program. 

The  calculation  of  the  potential  evapotranspiration  is  the  next 
operation  performed  by  the  main  program.   Potential  evapotranspiration  is 
determined  by  means  of  the  Penman  Combination  equation.   Following  comple- 
tion of  this  operation,  lake  and  bare  soil  evaporation  potentials  are 
determined  for  the  geographic  area  under  consideration.   Subroutine  SNOWRT 
is  then  called  by  the  main  program.   SNOWRT  evaluates  the  maximum  and 
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Figure  4.   Generalized  Flowchart  for  Subroutine  CROPCO 
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Figure  5.   Generalized  Flowchart  for  Subroutine  WTRMOD 
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minimum  daily  temperatures,  and  determines  if  the  precipitation  is  in  the 
form  of  rain  or  snow.   Also  the  potential  accumulation  of  the  snow  pack  or 
the  possibility  of  the  formation  of  a  snow  melt  is  evaluated.   Following 
these  computations  subroutine  SNOWRT,  represented  by  the  flowchart  shown 
in  Figure  6,  directs  the  computer  to  return  to  the  main  program. 

The  main  program  evaluates  the  amount  of  soil  moisture  within  the 
land  disposal  area  next.   If  any  precipitation  has  occurred,  the  magnitude 
of  infiltration,  interception  storage,  and  runoff  is  calculated.   The 
amount  of  moisture  which  infiltrated  from  the  surface  is  distributed 
within  the  upper  six  feet  of  the  soil  profile.   Any  moisture  infiltrating 
below  six  feet  is  considered  as  groundwater  recharge.   The  volume  for 
surface  runoff  is  calculated  by  using  the  method  developed  by  the  Soil 
Conservation  Service  (SCS)  for  the  antecedent  moisture  condition  II. 
Modifications  to  antecedent  moisture  conditions  I  cr  III  are  then  made 
and  evaluated.   The  soil  moisture  which  has  been  modified  is  then  redis- 
tributed.  The  main  program  calls  subroutine  DARCRT  where  redistribution 
is  accomplished  by  utilizing  a  modified  form  of  Darcy's  equation  for 
unsaturated  one-dimensional  flow.   The  equation  requires  the  use  of 
hydraulic  conductivity  and  soil-water  potential  parameters  which  should 
be  determined  from  laboratory  testing.   Subroutine  DARCRT,  as  shown 
conceptually  in  Figure  7,  then  directs  the  computer  to  return  to  the  main 
program.   If  soil  laboratory  testing  data  is  not  available,  moisture 
redistribution  can  be  established  through  an  evaluation  of  field  capacity, 
as  discussed  in  Chapter  II. 

The  next  operation  conducted  by  the  main  program  is  the  calculation 
of  the  magnitude  of  evapotranspiration  likely  to  occur  under  field 
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conditions.   This  calculation  is  dependent  upon  ambient  soil  moisture, 
vegetation,  and  climatic  conditions.   The  soil  moisture  distribution  is 
then  reevaluated.   After  the  determination  of  the  evapotranspiration  is 
complete  the  possibility  of  using  the  stored  water  for  irrigation  is 
considered.   By  adjusting  soil  moisture  values  while  decreasing  the  pond 
volume,  irrigation  can  be  simulated  mathematically  for  specific  disposal 
areas. 

If  feedlot  and/or  municipality  wastes  are  to  be  considered,  then 
subroutine  STORAG  is  called.   This  subroutine  generates  the  volume  of 
potential  feedlot  runoff  as  inflow  to  the  storage  facility  in  addition  to 
considering  any  variations  in  climatic  conditions.   Additional  pond  loading 
resulting  from  the  inflow  of  municipal  waste  disposal  can  also  be  evaluated. 
A  conceptual  flowchart  for  subroutine  STOPAG  is  shown  schematically  in 
Figure  3. 

Upon  returning  from  subroutine  STORAG  the  computer  calculates  the 
surface  evaporation  from  the  storage  pond.   The  required  volume  of  the. 
storage  pond  is  reevaluated  by  incorporating  all  the  inflows  to  the 
system.   If  the  combined  inflows  produce  a  volume  greater  than  the 
maximum  volume  available  for  storage,  then  the  amount  of  overflow  is 
calculated  and  a  message  is  issued  to  the  user  which  indicates  the  base 
dimensions  required  to  contain  all  the  inflow.   The  main  program  then 
updates  the  daily  and  monthly  accounts  for  the  storage  facility  and 
disposal  areas.   A  written  summary  is  provided  at  the  end  of  each  set 
of  yearly  calculations  prior  to  continuing  on  to  the  next  year.   The 
final  operation  of  the  computer  program  is  to  write  an  average  annual 
summary  which  can  be  used  for  a  prediction  analysis. 
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Figure  3.   Generalized  Flowchart  for  Subroutine  STORAG 
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Capabilities 

This  section  describes  the  capabilities  of  the  computer  program  to 
simulate  actual  field  conditions.   The  possible  existence  of  varying 
climatic,  soil  and  vegetation  conditions  as  well  as  inflow  and  outflow 
sources  requires  that  the  program  possess  multiple  capabilities. 

The  model  simulates  climatic  conditions  by  utilizing  previously 
recorded  precipitation  and  temperature  patterns  on  a  daily  basis  for  a 
specified  geographic  location.   The  required  climatological  data  is 
obtained  from  a  magnetic  computer  tape  provided  by  the  National  Weather 
Service  Climatic  Center  located  in  Asheville,  North  Carolina.   Missing 
data  on  the  tape  is  automatically  adjusted  in  the  program  by  assuming 
the  previous  day's  value.   In  addition,  the  program  can  be  used  to  predict 
the  possible  effects  of  weather  modification  based  on  historical  precipi- 
tation patterns.   The  period  and  location  to  be  evaluated  is  controlled 
by  input  parameters. 

The  major  components  used  in  the  program  are  the  inflows  from  the 
feedlot  and/or  municipal  wastes,  in  addition  to  the  size  of  the  storage 
and  disposal  facilities  to  be  used.   The  program  is  capable  of  evaluating 
the  results  from  up  to  nine  disposal  areas  under  the  following  three 
conditions:   (1)  a  watershed  with  runoff  as  inflow  to  the  storage 
facility,  (2)  a  watershed  with  runoff  as  inflow  to  the  storage  facility 
which  is  to  be  used  as  a  major  irrigation  source,  and  (3)  a  disposal  area 
requiring  only  supplemental  irrigation.   The  components  to  be  evaluated 
in  any  simulation  model  is  regulated  by  means  of  the  input  data. 

Variations  in  terms  of  area,  soil  type  and  crop  for  each  disposal 
area  can  also  be  incorporated.   The  twelve  soil  types  xv-hich  can  be 
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evaluated  by  this  program  are  presented  in  Appendix  B.   The  variation 
in  crops  which  can  occur  are:   wheat,  grain  sorghum,  corn,  soybeans, 
pasture,  alfalfa,  and  fallow.   An  annual  crop  rotation  using  periodic 
fallow  conditions  can  be  analyzed  with  crop  cover  occurring  on  either 
the  odd  or  even  years  of  the  simulation  period. 

The  hydrological  effects  on  each  area  are  evaluated  using  the  opera- 
tions defined  earlier  in  the  Chapter.   The  program  is  capable  of  evaluating 
the  redistribution  of  soil  moisture  by  using  a  modified  one-dimensional 
Darcy  equation  using  unsaturated  flow  conditions.   This  option  can  be 
used  only  when  the  unsaturated  flow  characteristics  of  the  soil  are  known. 

Irrigation  can  be  rapidly  evaluated  by  using  the  program.   The 
irrigation  application  rate,  in  units  of  inches  per  day  over  the  proposed 
areas,  can  be  varied  by  input  parameters.   In  addition,  the  application 
of  irrigation  can  be  optimized  by  specifying  that  the  soil  moisture  within 
the  top  twelve  inches  must  be  below  a  specified  percentage  of  field 
capacity  before  irrigation  is  implemented. 

The  disposal  facility  should  be  designed  to  provide  one  hundred 
percent  control  of  inflow  so  that  no  overflow  occurs.   Incorporated  into 
the  program  is  a  written  statement  to  the  user  indicating  the  required  base 
dimensions  to  provide  for  one  hundred  percent  control.   Statistical  hydro- 
logical  data  is  also  available  to  the  user  for  use  in  prediction  analysis. 
Since  the  cost  of  irrigation  is  increasing,  the  computer  program  developed 
previously  and  modified  during  this  research  study  provides  an  adequate 
and  efficient  means  to  design  and  evaluate  performance  of  systems  for 
utilizing  runoff  and  wastes  as  a  renewable  irrigation  resource. 
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Data  Input 

This  section  of  the  Chapter  describes  the  format  specified  for  the 
data  to  be  entered  into  the  computer  program.   The  section  is  divided  into 
seven  subsections  with  each  subsection  describing  the  required  format 
necessary  for  a  particular  set  of  data.   All  numeric  data  must  be  right-hand 
justified  and  a  complete  data  set  is  required  for  a  successful  mode], 
analysis.   Presented  in  Table  II  is  a  guide  for  data  input. 

Location  and  Mo de 1  Identification 

The  first  three  data  cards  are  written  in  alphanumeric  notation  for 
the  purpose  of  aiding  in  identification  of  a  particular  analysis.   The 
first  card  is  used  to  describe  the  geographic  location,  stored  as  NAME 
0?  CITY  AND  STATE.   This  description  begins  in  column  20  and  may  contain 
a  maximum  of  20  alphanumeric  characters.   The  second  card  is  used  to 
indicate  if  a  weather  modification  technique  is  to  be  implemented.   The 
term  WEATH  MODIF  MODEL,  which  begins  in  column  2  with  20  characters, 
indicates  that  a  weather  modification  scheme  was  used.   If  no  weather 
modification  technique  was  used  then  the  card  must  be  left  blank.   The 
third  data  card  is  used  to  identify  the  possible  incorporation  of  municipal 
wastes  to  the  storage  pond.   This  is  represented  by  the  term  INFLOW  TO 
STORAGE  POND,  beginning  in  column  2  with  24  alphanumeric  characters.   A 
blank  card  will  indicate  no  municipal  wastes  are  to  be  considered. 

Simulation  Model  Parameters 

The  data  cards  described  in  this  subsection  utilize  the  NAMELIST 
feature  of  the  Fortran  IV  language.   The  first  column  of  the  data  card 
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utilizing  this  feature  is  required  to  be  blank.   The  character  '&'  must 
appear  in  column  2  followed  immediately  by  the  NAMELIST  name,  with  no 
embedded  blanks.   The  NAMELIST  name  is  followed  by  a  blank  space.   All 
input  items  must  appear  in  the  form  of  a  variable  =  constant,  with  each 
item  separated  by  a  comma.   The  list  of  items  is  terminated  with  the 
characters  'SEND'  which  allows  for  several  data  cards  to  be  used  as  input 
under  one  NAMELIST. 

The  first  data  card  featuring  a  NAMELIST  option  is  termed  SEED.   The 
variables  included  are  MODEL  and  WPCNT.   MODEL  is  represented  by  the 
integer  l,2,3or4.   When  MODEL  is  equal  to  1,  there  will  be  no  modification 
of  the  precipitation  data.   When  MODEL  is  equal  to  2,  the  'B-Zero'  modifi- 
cation, previously  described,  is  implemented.   When  3  is  used,  rainfall  is 
increased  by  the  percentage  inputted  for  WPCNT.   The  integer  4  indicates 
the  implementation  of  the  '3-Zero  Modified'  alternation  scheme. 

The  next  set  of  data  is  read  from  the  NAMELIST  ALPHA.   The  following 
variables  are  included  in  this  NAMELIST  feature:   BRUNTA,  BRUNTB,  E,  RCROP, 
DORM,  GROW,  PAVLU,  PCVMAX  and  DSRATE.   The  first  four  variables  are  used 
in  computing  the  evapotranspiration  potential.   BRUNTA  and  BRUNTB  are 
geographic  constants,  E  is  the  wind  coefficient,  and  RCROP  describes  the 
reflectance  coefficient  for  use  in  the  Penman  equation  described  previously. 
The  variables  DORM  and  GROW,  expressed  in  inches,  are  used  in  computing 
surface  funoff  from  a  feedlot  area. 

Scheduling  irrigation  requires  the  input  of  PAVLU,  PCVMAX,  and  DSRATE. 
Irrigation  is  implemented  only  when  the  soil  moisture  is  below  a  percentage, 
PAVLU,  of  field  capacity.   PCVMAX  expresses  the  minimum  pond  volume,  as 
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a  fraction  of  the  maximum  volume,  required  for  irrigation.   The  disposal 
rate,  in  units  of  inches  per  day,  is  defined  by  DSRATE. 

The  final  NAMELIST  data  set  is  represented  by  the  name  BETA.   The 
variables,  HMAX,  L,  W,  and  S  are  the  base  parameters,  in  units  of  feet, 
for  maximum  depth,  length,  width  and  side  slope  used  in  sizing  the 
retention  pond.   Variables  INDST,  YSTART,  MSTART,  and  YEND  locate  the 
appropriate  block  of  precipitation  and  temperature  data  to  be  read  from 
the  magnetic  tape.   INDST  is  an  index  reference  number  for  the  geographic 
location.  YSTART,  MSTART  and  YEND  define  the  starting  year,  month,  and 
ending  year  of  the  simulation  period. 

The  variable  ROTAYR  is  represented  by  the  integer  1  or  2.   When 
ROTAYR  is  equal  to  1,  an  annual  crop/fallow  rotation  is  implemented  with 
the  crop  analysis  on  the  first  year  of  the  simulation  period.   When  ROTAYR 
is  equal  to  2,  fallow  is  established  on  the  first  year  of  the  simulation 
period. 

DARCEQ  is  also  represented  by  the  integer  1  or  2.   When  DARCEQ  is 
equal  to  1,  the  redistribution  of  moisture  within  the  soil  profile  is 
evaluated  by  the  Darcy  equation  for  one-dimensional  unsaturated  flow. 
When  DARCEQ  is  equal  to  2,  the  simplified  method,  described  previously, 
is  used  for  analyzing  the  redistribution  of  moisture.   The  unsaturated  flow 
data  cards,  described  below,  are  not  required  when  DARCEQ  is  equal  to  2. 

The  variable  BYPASS  establishes  feedlot  runoff  and/or  municipal  waste 
as  a  source  of  inflow  to  the  storage  facility.   BYPASS  is  represented  by 
the  integer  l,2,3or4  with  each  integer  establishing  the  following  analysis: 
(1)  BYPASS  equals  1-  no  feedlot  runoff  or  municipal  waste,  (2)  BYPASS 
equals  2-  feedlot  runoff,  (3)  BYPASS  equals  3-  municipal  waste,  and 
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(4)  BYPASS  equals  4-  both  feedlot  runoff  and  municipal  waste.  Data  cards 
for  municipal  waste  input,  described  below,  are  required  only  when  BYPASS 
is  greater  than  or  equal  to  3. 

The  final  three  parameters  included  in  BETA  are  LTAREA,  STORM,  and 
NPLOTS.   The  feedlot  area,  in  acres,  is  represented  by  LTAREA.   The 
statistical  25  year-24  hour  storm  is  stored  as  STORM,  in  inches.   NPLOTS 
indicates  the  number  of  disposal  areas  included  in  the  analysis. 


Meteorological  Data 


The  data  required  for  the  meteorological  conditions  is  entered  on  the 
next  set  of  data  cards.   This  data  is  entered  on  a  monthly  basis  with 
each  month  appearing  on  a  separate  card,  thus  a  total  of  twelve  cards  are 
necessary  to  describe  the  expected  conditions. 

The  mean  monthly  sunshine  and  relative  humidity  are  entered  in 
columns  3  and  5,  respectfully.   Both  variables  are  entered  in  units  of 
percent  and  are  stored  by  the  variables  PSUNS  and  RHD.   RA,  the  mid- 
monthly  intensity  of  solar  radiation  in  mm  of  water  evaporated  per  day, 
appears  in  column  7.   In  column  11  the  mean  monthly  wind  speed,  WIND, 
is  entered  in  the  units  of  miles  per  hour.   The  final  parameter  is  stored 
as  MMAT,  the  mean  monthly  temperature  in  degrees  Farenheit,  and  is  located 
in  column  14. 

Disposal  Area  Data 

The  next  group  of  data  requires  one  card  for  each  disposal  area  under 
analysis.   The  first  card  is  read  in  matrix  form  (i,j)  where  i  and  j 
represent  the  reference  number  for  crop  and  soil  in  column  3  and  6, 
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respectfully.   This  information  is  stored  by  the  computer  under  the 
parameter  name  IAREA.   The  area  of  each  plot,  stored  as  AREA,    is  entered 
on  the  second  card  in  column  2,  using  units  of  acres.   The  parameter 
IPLAN  on  the  third  card  indicates  the  reference  number  for  irrigation 
or  runoff  collection  for  each  disposal  area.   IPLAN  is  an  integer 
located  in  column  2.   Column  3  of  this  same  card  implements  annual  crop 
rotations  through  an  integer  reference  number  stored  as  ROTATE.   All 
reference  numbers  are  defined  in  Appendix  B. 

The  final  data  card  establishes  the  dates  for  planting  and  harvesting 
crops  on  each  area.   MGSBP  and  DGSBP  indicate  the  month  and  date  for  the 
beginning  of  the  growing  season,  located  in  column  2  and  5,  respectfully. 
The  ending  month  and  date  are  stored  as  MGSEP  and  DGSE?  in  columns  8  and  11, 
respectfully.   The  monthly  parameters  are  entered  as  integers  with  January 
represented  by  1,  February  by  2,  etc.   The  field  for  each  integer  is  two 
spaces. 

Unsaturated  Flow  Data 

A  data  card  for  the  hydraulic  properties  of  each  disposal  area  is 
required  when  the  moisture  redistribution  is  to  be  analyzed  by  the  Darcy 
equation  for  one-dimensional  unsaturated  flow.   The  parameters  CONDUZ, 
C0NDL2,  CONDGW,  entered  in  columns  1,  10,  and  20,  respectfully,  are  the 
computed  hydraulic  conductivities  for  the  soil  profile  layers  described 
previously.   The  soil-water  potential  variables,  stored  as  HUZ,  HLZ,  HGWZ, 
are  entered  in  columns  30,  40,  and  50,  respectfully. 
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Municipal  Waste  Data 

The  municipal  waste  data  is  entered  on  a  monthly  basis,  requiring 
a  total  of  12  cards.   The  monthly  values,  stored  as  MUNDIS,  are  located 
in  column  1.   MUNDIS  is  in  the  units  of  gallons  per  day  with  each  value 
being  the  average  monthly  disposal.   These  data  cards  are  included  only 
when  municipal  waste  is  included  as  input  to  the  storage  facility. 

Climatological  Data 

The  final  set  of  data  cards  are  IBM  control  cards  generating  the 
input  of  climatic  data  from  a  magnetic  tape.   The  tapes,  provided  by  the 
National  Weather  Service  Climatic  Center,  are  written  in  NWS  Format  II 
language.   The  data  includes  daily  precipitation,  inches,  and  temperature, 
degrees  Farenheit,  for  the  simulation  period. 

Data  Output 

This  section  of  the  Chapter  describes  the  output  data  from  the 
computer  program.   The  section  is  divided  into  three  subsections  with 
each  describing  a  portion  of  the  program  output.   An  example  of  the  output 
data  is  presented  in  Appendix  E. 

Input  Data 

The  first  set  of  data  to  appear  as  output  is  a  summary/  of  input 
parameters.   The  first  line,  indicating  the  station  and  simulation 
period,  is  followed  by  the  model  identifiers  and  the  size  of  the  25  year- 
24  hour  storm.   The  retention  pond  variables  are  then  listed.   Included 
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are  the  facilities  base  dimensions,  side  slope,  maximum  depth,  maximum 
pond  volume  and  the  receiving  area  for  precipitation. 

Next  the  area  variables  are  presented  according  to  plot  numbers. 
Under  each  plot  the  area,  soil  type  and  vegetative  cover,  including  an 
annual  crop  rotation,  is  listed.   The  plan  to  be  implemented  on  each  area 
(irrigation  and/or  runoff  collection)  and  the  irrigation  management  scheme 
is  also  included. 

The  summary  of  input  parameters  described  under  this  subsection  is 
repeated  as  output  prior  to  the  final  summary.   The  purpose  for  duplicating 
this  information  is  to  aid  the  user  in  analyzing  the  results. 

Annual  Summary 

This  subsection  describes  the  annual  summaries  for  the  disposal 
facility  and  areas  during  each  year  of  the  simulation  period.   A  message 
to  the  user  is  initiated  at  the  beginning  of  the  annual  summary  when  the 
following  events  occur:   (1)  a  storm  of  magnitude  equal  to  or  greater  than 
the  25  year-24  hour  storm,  and  (2)  an  overflow  of  the  storage  facility. 
The  date  of  occurrence  and  also  the  dimensions  required  to  prevent  overflow 
are  included  in  the  statements. 

The  water  budget  for  the  storage  facility  is  printed  in  tabular 
form.   Each  row  tabulates  the  monthly  summaries  with  the  last  row  repre- 
senting the  annual  totals.   The  units  incorporated  for  inflow  and  outflow 
are  acre- inches.   Columns  2  through  5  correspond  to  inflows  to  the  system 
while  columns  7  through  9  represent  outflows.   The  sources  of  inflow 
include  precipitation,  municipal  waste,  feediot  runoff  and  watershed 
runoff.   Irrigation, surface  evaporation  and  overflow  are  the  sources  of 
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outflow.   In  addition,  the  number  of  days  irrigation  is  implemented  and 
the  depth,  in  feet,  of  the  storage  facility  at  the  end  of  each  month  is 
included  in  volumns  6  and  11,  respectfully.   Column  10  represents  the  change 
in  volume  in  units  of  acre- inch.   This  column  can  be  used  as  a  water  account 
check  since  the  inflows  minus  outflows  equals  the  change  in  volume. 

Next  the  moisture  budget  for  each  area  is  printed  in  tabular  form. 
Sources  of  moisture  input  for  each  area  include  precipitation  and  irrigation, 
The  outputs  are  interception,  surface  runoff,  percolation  and  evapotrans- 
piration.   A  column  indicating  the  total  monthly  and  annual  change  in  soil 
moisture  is  also  included.   All  units  are  in  terms  of  inches.   The  moisture 
account  balances  when  inputs  minus  outputs  minus  change  in  snow  storage 
equals  change  in  soil  moisture. 

Following  the  monthly  accounts  for  each  area  an  annual  summary  of  the 
operating  parameters  is  printed.   This  includes  statements  indicating  the 
percent  of  wastewater  controlled  and  percent  of  maximum  pond  volume  required, 
Also  included  is  the  accumulation  of  a  snow  pack  during  the  season  and  the 
change  in  snow  storage  as  compared  to  the  previous  year. 

Final  Summary 

The  final  portion  of  output  is  a  final  summary  utilizing  average 
annual  computations  over  the  simulation  period.   The  summaries  include 
meteorological  input,  pond  operations  and  individual  disposal  area 
averages  which  can  be  used  in  a  cost  to  benefit  prediction  analysis. 
Also  included  is  a  summary  of  statistical  data  used  to  generate  the 
precipitation  and  runoff  frequency  curves. 


CHAPTER  V 

RESULTS  AND  DISCUSSION 

This  Chapter  presents  and  discusses  the  results  obtained  from  using 
the  computer  program  to  analyze  the  effects  of  weather  modification  and 
the  design  of  a  supplemental  irrigation  system  using  municipal  wastewater. 
The  first  analysis  presented  is  the  evaluation  of  the  hydrological  effects 
of  implementing  a  weather  modification  program  in  Kansas.   The  second 
analysis  consists  of  an  optimized  design  for  using  effluent  municipal 
wastewater  for  supplemental  irrigation. 

Weather  Modification 

Ten  locations  representing  a  geographic  cross-section  of  Kansas  were 
used  to  test  the  effects  of  cloud  seeding  on  the  current  hydrological 
environment  within  the  state.   The  locations  selected  to  represent  the 
western  section  of  Kansas  were  Colby,  Dodge  City,  Garden  City,  Goodland 
and  Hays.   Belleville  and  Ellsworth  were  selected  to  represent  the  central 
section  of  the  state.   The  eastern  portion  was  represented  by  Horton, 
Independence  and  Topeka.   These  ten  locations  were  chosen  because  of  the 
availability  of  historical  climatic  data  in  these  regions.   Each  location 
was  analyzed  for  a  crop  cover  of  wheat  without  irrigation.   The  period 
simulated  for  each  location  used  in  the  analysis  is  presented  in  Table  III. 
The  input  data  used  for  each  location  is  presented  in  Appendix  D. 
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TABLE  III 
SUMMARY  OF  RESULTS  OF  PRECIPITATION  AND  RUNOFF 


Prec 

ipitation 

(in.) 

Runoff  (in 

.) 

Time 

B-ZERO 

B-ZERO 

Location 

Period 

Natural 

B-ZERO 

MODIFIED 

Natural 

B-ZERO 

MODIFIED 

Colby 

1950-1962 

18.4 

20.9 

21.3 

0.8 

0.9 

1.1 

Dodge  City 

1949-1973 

20.7 

23.0 

23.6 

1.1 

1.2 

1.4 

Garden  City 

1950-1974 

18.2 

20.2 

20.7 

0.3 

0.8 

1.0 

Goodland 

1949-1973 

15.9 

13.4 

18.7 

0.4 

0.5 

0.6 

Hays 

1948-1973 

23.4 

25.4 

26.2 

1.6 

1.6 

2.0 

Belleville 

1949-1973 

30.0 

32.1 

33.2 

3.4 

3.5 

4.0 

Ellsworth 

1946-1970 

27.1 

29.1 

30.0 

2.4 

2.5 

2.9 

Ho r ton 

1946-1970 

36.7 

33.6 

40.1 

5.9 

5.9 

6.8 

Independence 

1948-1961 

36.9 

38.8 

40.0 

5.7 

5.7 

6.4 

Topeka 

1964-1973 

37.0 

33.7 

40.1 

5.6 

....  .  . 

5.5 

6.4 
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Precipitation 

Two  weather  modification  schemes  were  used  to  modify  the  historical 
precipitation  patterns  for  each  of  the  ten  locations  during  the  months 
March  through  September.   The  first  modifying  scheme,  Model  B-ZERO,  produces 
a  seventy-five  percent  increase  in  precipitation  from  clouds  producing 
normally  0.10  inch  or  less.   Precipitational  clouds  which  normally  pro- 
vide between  0.11  and  0.50  inches  are  increased  thirty  percent  under 
Model  B-ZERO  conditions  while  amounts  between  0.51  and  1.00  inches  are 
increased  ten  percent.   Large  precipitational  clouds  are  reduced  ten 
percent  by  Model  B-ZERO. 

The  second  weather  modification  scheme  analyzed  is  Model  B-ZERO 
MODIFIED.   This  model  modifies  the  daily  rainfalls  in  Kansas  the  indicated 
percentages  given  for  Model  B-ZERO  with  the  exception  of  large  rainfalls. 
There  is  no  modification  to  rainfalls  of  greater  than  1.00  inch  under 
this  model.   Model  B-ZERO  MODIFIED  provides  the  potential  range  in 
precipitation  increase  using  cloud  seeding  programs. 

The  results  of  modifying  the  daily  rainfalls  during  the  periods 
simulated  for  each  location  are  presented  in  Table  III.   The  average 
increase  in  the  magnitude  of  precipitation  for  the  time  period  considered 
within  the  state  was  approximately  two  inches  per  year  using  Model  3-ZER0 
while  approximately  three  inches  per  year  was  obtained  using  Model  B-ZERO 
MODIFIED.   These  results  agree  with  previous  research  estimates  which 
indicated  that  a  possible  increase  of  ten  percent  could  be  obtained  from 
using  a  cloud  seeding  program  (7) . 

Although  the  increase  in  magnitude  of  rainfall  was  approximately 
equal  throughout  the  state,  the  percentage  of  rainfall  increase  varied 
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significantly  for  western,  central  and  eastern  portions  of  Kansas. 
In  the  western  section  of  the  state,  Model  B-ZERO  indicated  an  average 
increase  in  precipitation  of  11.6  percent  while  B-ZERO  MODIFIED  indicated 
a  14.5  percent  increase.   An  increase  of  7.3  and  10.7  percent  resulted  in 
the  central  portion  of  Kansas  from  using  Models  B-ZERO  and  B-ZERO  MODIFIED 
respectfully.   Increases  in  the  eastern  section  were  indicated  to  be  only 
5.0  percent  for  B-ZERO  and  8.7  percent  for  B-ZERO  MODIFIED.   The  range  of 
this  potential  precipitational  increase  is  shown  schematically  in  Figures 
9,  10  and  11  for  the  western,  central  and  eastern  sections  respectfully. 
The  time  period  chosen  was  from  1950  through  1962  because  of  the  wealth 
of  data  available  during  this  period  for  the  ten  locations. 

Surface  Runoff 

The  surface  runoff  which  could  be  expected  was  estimated  using  the 
method  developed  by  the  Soil  Conservation  Service  previously  described 
in  Chapter  II.   It  should  be  noted  that  prior  to  any  surface  runoff  an 
interception-storage  requirement  of  0.1  inch  had  to  be  satisfied.   The 
results  obtained  for  the  surface  runoff  based  on  the  ten  locations 
previously  defined  are  presented  in  Table  III. 

As  can  be  seen  in  Table  III  the  expected  average  amount  of  annual 
surface  runoff  in  western  Kansas  is  only  0.9  inches,  while  central  and 
eastern  Kansas  averaged  2.9  and  5.7  inches  respectfully.   The  use  of 
Model  B-ZERO  increased  the  runoff  in  all  three  sections  of  Kansas  approxi- 
mately 0.1  inch.   However,  Model  B-ZERO  MODIFIED  indicated  an  increase  in 
the  annual  runoff  of  0.2  inches  in  western  Kansas,  0.5  inches  in  central 
Kansas  and  0.3  inches  in  eastern  Kansas.   It  is  believed  that  the  results 
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Figure  9.   Natural  and  Enhanced  Annual  Precipitation 
for  Western  Kansas 
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Figure  10.   Natural  and  Enhanced  Annual  Precipitation 
for  Central  Kansas 
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Figure  11.   Natural  and  Enhanced  Annual  Precipitation 
for  Eastern  Kansas 
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obtained  from  using  the  two  precipitation  enhancement  schemes  accurately 
describe  the  expected  range  of  increase  in  runoff.   This  belief  is  based 
on  the  results  of  the  research  conducted  by  Hanson  and  Woolhiser  (18) . 
Their  research  conducted  on  conditions  in  western  South  Dakota  which  were 
similar  to  those  used  in  this  study  indicated  a  limited  increase  in 
runoff  similar  to  3-ZERO  and  a  27  percent  increase  using  a  Model  similar 
to  B-ZERO  MODIFIED. 

Deep  Percolation 

The  deep  percolation  values  presented  in  Table  IV  represent  the 
quantity  of  moisture  passing  through  the  lower  zone  of  the  soil  profile 
as  discussed  previously.   The  natural  amount  of  deep  percolation  under 
normal  precipitational  patterns  in  the  western  section  of  the  state  is 
apparently  quite  low.   The  computer  analysis  predicts  an  increase  in  the 
magnitude  of  deep  percolation  to  be  0.1  inch  annually  using  a  precipita- 
tion enhancement  program  in  western  Kansas.   The  central  section  of  the 
state,  which  currently  has  an  average  annual  deep  percolation  of  0.7 
inches  could  expect  an  increase  of  0.5  to  0.8  inches. 

The  eastern  portion  of  Kansas  would  receive  the  greatest  benefit  to 
potential  groundwater  recharge  as  a  result  of  using  a  weather  modification 
scheme.   The  average  amount  of  deep  percolation  under  natural  conditions 
is  currently  estimated  to  be  3  inches  annually.   Models  B-ZERO  and  B-ZERO 
MODIFIED  predict  that  potential  recharge  could  be  increased  to  4.0  and  4.5 
inches  respectfully.   These  values  indicate  that  approximately  50  percent 
of  any  additional  moisture  which  would  result  from  cloud  seeding  in  eastern 
Kansas  could  be  considered  as  potential  groundwater  recharge. 
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TABLE  IV 
SUMMARY  OF  RESULTS  FOR  RECHARGE  AND  EVAPOTRANSPIRATION 


Location 


Time 
Period 


Deep  Percolation  (in.) 


B-ZERO 


Evapotranspiration  (in.) 


B-ZERO 


Natural  3-ZERO  MODIFIED  Natural   3-ZERO  MODIFIED 


Colby 

Dodge  City- 
Garden  City 
Goodland 
Hays 

Belleville 
Ellsworth 
Horton 
Independence 
Topeka 


1950-1962 
1949-1973 
1950-1974 
1949-1973 
1948-1973 
1949-1973 
1946-1970 
1946-1970 
1948-1961 
1964-1973 


0 

0 

0 

0 
0.10 
1.24 
0.22 
4.25 
2.06 
3.01 


0 
0.04 
0.03 

0 
0.14 
1.99 
0.46 
5.51 
2.95 
3.63 


0 
0.07 
0.05 

0 
0.18 
2.36 
0.61 
6.05 
3.33 
4.05 


12.9 
14.5 
12.9 
10.8 
16.6 
19.3 
18.9 
19.8 
22.1 
21.2 


14.8 
16.2 
14.4 
12.7 
18.2 
20.1 
20.2 
20.2 
22.7 
21.9 


15.0 
16.6 
14.3 
12.9 
13.6 
20.3 
20.6 
20.2 
22.9 
22.0 
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Actual  Evapo t ransp irat ion 

Evapotranspiration  (AET)  as  used  in  this  study  includes  both  bare  soil 
evaporation  and  any  transpiration.   The  AET  values  predicted  for  each  loca- 
tion are  presented  in  Table  IV.   Using  the  data  presented  in  this  Table  the 
average  annual  AET  values  were  estimated  to  be  13.5,  19.1  and  21.0  inches 
for  western,  central  and  eastern  Kansas  respectfully.   These  values  indi- 
cate that  nearly  two-thirds  of  the  total  precipitation  under  natural  con- 
ditions is  dissipated  by  AET. 

In  western  Kansas,  Model  B-ZERO  predicted  an  increase  in  the  AET 
value  of  1.8  inches  while  Model  B-ZERO  MODIFIED  indicated  an  increase 
of  2  inches.   These  values  imply  that  75  percent  of  the  moisture  added 
as  a  result  of  weather  modification  will  be  utilized  in  the  ET  process. 
In  the  central  portion  of  the  state,  increases  of  1.0  to  1.3  inches  were 
predicted.   These  values  represent  nearly  50  percent  of  the  precipitation 
added  from  a  cloud  seeding  program.   The  increased  values  in  evapotrans- 
piration suggest  an  economic  benefit  through  potential  improvement  of 
crop  yields  due  to  weather  modification  in  western  and  central  Kansas. 
Eastern  Kansas  is  predicted  to  have  a  maximum  AET  increase  of  0.7  inches 
which  suggests  that  the  soil  moisture  required  for  maximizing  evapotrans- 
piration is  adequate  in  this  portion  of  the  state. 

Municipal  Wastewater  Effluents 

This  section  of  the  Chapter  presents  the  results  obtained  by  the 
computer  program  to  analyze  and  design  a  supplemental  irrigation  system 
using  the  effluent  from  a  municipal  wastewater-treatment  plant.   The 
application  of  municipal  wastewater  as  a  supplement  to  an  irrigation 
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system  requires  the  design  of  a  storage  facility.   Since  the  volume  of 
municipal  wastewater  is  generally  large  and  the  flow  rate  is  time 
dependent  the  storage  facility  must  act;  as  a  fluid  control  system  in 
addition  to  being  based  on  an  optimum  design. 

For  comparative  purposes  two  designs  were  made  for  a  storage 
control  facility  to  be  located  in  the  Topeka  area  during  a  ten  year 
simulation  period  beginning  in  1964  and  ending  in  1973.   The  first  design 
was  for  a  storage  facility  having  the  following  natural  sources  of 
inflow:   (1)  precipitation  directly  into  the  facility,  and  (2)  surface 
runoff  from  a  120  acre  watershed  area.   The  second  design  included 
municipal  inflow  into  the  storage  facility  using  an  average  flow  race 
of  84,000  gallons  per  day.   Both  facilities  were  anaylized  for  supplying 
irrigation  water  for  two  40  acre  areas  with  the  vegetative  cover  of  corn 
and  wheat.   Results  are  presented  in  Table  V. 

The  use  of  municipal  wastes  as  supplemental  irrigation  water  required 
a  storage  facility  volume  of  2720  acre-inches  for  controlling  ninety-eight 
percent  of  the  inflow  while  under  the  natural  inflow  conditions  the  storage 
volume  required  was  690  acre- inches.   Utilizing  the  municipal  effluent 
the  total  annual  increase  in  the  volume  of  irrigation  water  applied  to 
the  disposal  areas  was  838  acre- inches.   Irrigation  was  applied  at  a 
disposal  rate  of  0.5  inches  per  day  when  the  soil  moisture  was  below  ninety 
percent  of  field  capacity. 

The  additional  irrigation  increased  potential  groundwater  recharge 
an  average  of  4.6  inches  annually.  Evapotranspiration  on  the  irrigated 
corn  field  increased  2.9  inches  annually.   On  the  wheat  field 


TABLE  V 
SUMMARY  OF  RESULTS  FOR  A  DISPOSAL  CONTROL  FACILITY 
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Parameter 


Natural  Inflow 


Natural  and 
Municipal  Inflow 


Storage  Volume 

690  ac-in 

2720  ac-in 

Receiving  Area 

7.1  ac 

26.5  ac 

Precipitation 

38.7  in 

38.7  in 

Moisture  Deficit 

3.5  in 

3.5  in 

Disposal  Area  -  Corn 

Recharge 

5.1  in 

9.4  in 

No.  Irrigation  Days 

17 

40 

Evapo transpiration 

25.0  in 

27.9  in 

Disposal  Area  -  Wheat 

Recharge 

6.7  in 

11.6  in 

No.  Irrigation  Days 

14 

33 

Evapo transpiration 

22.7  in 

23.8  in 
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evapotranspiration  was  increased  1.1  inches  annually.   The  increase  in 
evapotranspiration  indicates  the  possibility  of  maximizing  crop  production, 

The  computer  program  provides  an  effective  means  of  evaluating 
the  application  of  municipal  wastewater  effluent  on  agricultural  land. 
To  optimize  the  storage  facility  design  various  disposal  rates  as  well 
as  differing  soil  conditions  should  be  considered. 


CHAPTER  VI 
CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

The  objective  of  this  research  study  was  to  develop  a  new  methodo- 
logical approach  to  the  problem  of  providing  a  dependable  supply  of  water 
for  man's  use.   This  objective  was  accomplished  by  developing  and  testing 
a  multipurpose  hydrologic  yield  model.   The  model  was  used  to  evaluate 
the  effects  of  weather  modification  schemes  on  the  hydrological  environ- 
ment in  Kansas.   In  addition,  the  utilization  of  municipal  waste  as  a 
controllable  source  of  water  for  irrigation  was  evaluated. 

As  a  result  of  this  study,  the  following  conclusions  can  be  made: 

1.  Cloud  seeding  programs  in  Kansas  could  produce  a  uniform  increase 
in  the  total  amount  of  precipitation. 

2.  Weather  modification  could  be  used  to  increase  surface  runoff, 
deep  percolation  and  evapo transpiration. 

3.  Fifty  percent  of  any  moisture  added  as  a  result  of  cloud  seeding 
in  eastern  Kansas  could  be  used  to  recharge  the  groundwater  supplies. 

4.  In  western  and  central  Kansas  more  than  fifty  percent  of  any 
moisture  added  because  of  cloud  seeding  could  be  used  in  the  evapotrans- 
piration  process. 

5.  The  computer  program  provides  an  effective  method  of  evaluating 
the  hydrological  effects  which  result  from  the  implementation  of  a  weather 
modification  scheme. 
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6.   The  computer  program  provides  an  effective  means  of  predicting 
the  optimized  design  for  a  disposal  control  facility. 

Recommendat  ions 

With  respect  to  future  research,  the  following  recommendations  are 
made: 

1.  Additional  verification  of  the  simulation  model  should  be 
conducted  under  both  controlled  laboratory  conditions  and  field  conditions. 

2.  A  comprehensive  study  should  be  conducted  to  evaluate  moisture 
redistribution  utilizing  Darcy's  equation  for  unsaturated  one-dimensional 
flow. 

3.  A  comprehensive  study  should  be  undertaken  to  evaluate  municipal 
wastewater  as  a  controllable  source  of  irrigation  water  for  various 
geographic  and  climatic  conditions. 

4.  A  comprehensive  study  should  be  initiated  to  study  the  possibility 
of  collecting  surface  runoff  from  natural  watersheds  as  a  supplemental 
irrigation  water  source. 

5.  The  computer  program  should  be  modified  to  incorporate  the 
following:   (1)  seepage  from  the  storage  facility,  (2)  the  possibility 
of  constructing  multiple  storage  facilities,  and  (3)  a  cost  optimization 
package. 
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AET  actual  evapotranspiration 

E  soil  evaporation 

ES  saturated  vapor  pressure 

ETP  evapotranspiration  potential 

K  hydraulic  conductivity 

N  runoff  curve  number 

P  precipitation 

PSUNS  percent  sunshine 

Q  surface  runoff 

Ra  solar  radiation 

RHD  relative  humidity 

S  maximum  potential  difference 

T  mean  daily  air  temperature 

Ta  mean  daily  air  temperature 

WVD  wind  speed 

a  geographic  location  constant 

b  geographic  location  constant 

c1  hydraulic  coefficient 

e  mass  transfer  coefficient 

i  hydraulic  gradient 

k  crop  consumptive  use  coefficient 

q?  volume  of  water  movement  per  unit  area 

r  reflectance  coefficient 

t  time  after  first  stage  evaporation 
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Ah  change  in  soil  water  potential 

At  time  increment 

Az  change  in  elevation 

0„  available  soil  moisture 
a 

G  maximum  available  soil  moisture 

max 
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TABLE  VI 
RUNOFF  CURVE  NUMBERS  FOR  ANTECEDENT  SOIL  MOISTURE  CONDITION  II* 


Crop 

Irrigation  ■ 

Soil 

Grain 

Fallow 

Class 

Wheat 

Sorghum 

Corn 

Soybeans 

Pasture 

Alfalfa 

Soil 

1 

33 

86 

36 

84 

80 

86 

34 

2 

83 

86 

86 

84 

80 

86 

84 

3 

78 

82 

82 

81 

74 

82 

81 

4 

78 

82 

82 

81 

74 

82 

81 

5 

69 

75 

75 

73 

61 

75 

73 

6 

69 

75 

75 

73 

61 

75 

73 

7 

69 

75 

75 

73 

61 

75 

73 

8 

69 

75 

75 

73 

61 

75 

73 

9 

69 

75 

75 

73 

61 

75 

73 

10 

69 

75 

75 

73 

61 

75 

73 

11 

69 

75 

75 

73 

61 

75 

73 

12 

55 

65 

65 

61 

39 

65 

61 

*From  Bean,  1976 
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TABLE  VII 

SOIL  CLASSIFICATION  USED  FOR  IRRIGATION  DESIGN* 


Profile 

Irrigation   Depth 

Soil  Class    (ft.)  Soil  Class  Description 

1         3'        Deep  soils  with  silt  loam  or  silty  clay  loam 

surface  layers  and  slowly  to  very  slowly  permeable 
heavy  clay  and  claypan  subsoils. 

3'        Deep  soils  with  silty  clay  or  clay  textures 

throughout.   Surface  infiltration  and  subsoil 
permeability  are  very  slow  when  the  soil  is 
moist.   Shrinkage  from  drying  causes  extensive 
cracking,  resulting  in  high  infiltration  rates 
until  swelling  occurs. 

3  5'        Deep  soils  with  silt  loam,  loam,  clay  loam,  or 

silty  clay  loam  surface  layers  and  clay  loam, 
silty  clay  loam,  or  silty  clay  subsoils. 
Subsoil  permeability  is  slow  to  moderately  slow. 
Shrinkage  cracks  resulting  from  drying  in  the 
soils  with  more  clayey  subsoil  textures  give  a 
relatively  high  initial  infiltration  rate. 

4  2.5'        Moderately  deep  soils  with  silt  loam,  clay  loam, 

or  silty  clay  loam  surface  layers  and  clay  loam 
or  silty  clay  subsoils  with  predominately 
moderately  slow  permeability. 

5  5*        Deep  soils  with  silt  loam, loam,  clay  loam,  or 

silty  clay  loam  surface  layers  and  subsoils. 
Subsoil  permeability:   moderate  to  moderately 
slow. 

6  3'        Moderately  deep  soils  with  silt  loam  or  loam 

surface  layers  and  loam,  clay  loam,  or  silty 
clay  loam  subsoils  with  moderate  to  moderately 
slow  permeability. 

7  5'        Deep  soils  with  silt  loam,  loam  or  very  fine 

sandy  loam  surface  layers  and  moderately  permeable, 
medium  textures  subsoils. 

8  2.5'        Moderately  deep  soils  with  silt  loam,  loam  or 

very  fine  sandy  loam  surface  layers  and  moderately 
permeable  clay  loam,  loam,  or  silt  loam  subsoils. 
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9         5'        Deep  soils  with  fine  sandy  loam  and  loam  surface 

layers  and  subsoils  that  have  moderately  rapid 
permeability.   Available  water  capacity  is 
moderate  to  low. 

10  5'        Soils  are  moderately  deep  over  sand  with  sandy 

loam  to  loam  surface  layers  and  moderately  rapid 
to  rapidly  permeable  subsoils  with  low  available 
water  capacity. 

11  5'        Deep  soils  with  loam  fine  sand  or  loamy  sand 

surface  layers  and  moderately  rapid  to  rapidly 
permeable  subsoils. 

12  5'        Deep  rapidly  permeable  soils  with  sand  or  fine 

sand  textures  throughout. 


*From  Bean,  1976 
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TABLE  VIII 
CROP  CODES  USED  IN  THE  COMPUTER  PROGRAM 


Crop 


Reference  Number 


Wheat 

1 

Grain  Sorghum 

2 

Corn 

3 

Soybean 

4 

Pasture 

5 

Alfalfa 

6 

Fallow  Soil 

7 

TABLE  IX 
AREA  CODES  USED  IN  THE  COMPUTER  PROGRAM 


Plan 
Indicated 


Rotation 
Indicated 


Reference  Number 


Collect  Area  Runoff 

Collect  Area  Runoff 
and  Irrigate 

Irrigate 


Crop  -  Fallow 
None 
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FORTRAN    IV    G    LEVEL      21  MAIN  OATE    «     79045  21/27/40 

C 

C  HYORGLOGICAL  HATERSHEO  MOCEL 

C 

C  KANSAS  STATE  UNIVERSITY    197  8 

C 

0001  INTEGER    CRCP. FROZE , SO  I L ,STI NJ  ,T , Y EAR , YE ARS , YENC.Y ST  ART 

0002  INTEGER    PL CT, PREVYK , 3YPAS S, NR OPTS , DAKCEQ ,KUTA YR 

0003  INTEGER    DGS8, DCS t , DCS GP(9  ) , CGSEP I  9 ) , C AY , RCT AT E I  9 ) 

0034  REALIA,  I  AAOU,  IAFT,KCRQP,KS,LAKEVP,L.<EVPI,LTAREA,L.M,MA,Nn,NOPERC 

0C05  HEAL    LI  , MM  AT  I  12),  IRRVCL,  I  RRSUM.MAXVOL  ,MR,NRNOF  ,LZSM,MUNOI  S(  12) 

0006  DIMENSION    HUZ(9) ,HLZ I  9) ,HCw < 9 ) , EC < 9 ) ,  T  (  9 )  .SMPCI9) .DISVGL19) 

000  7  DIMENSION    SM (9  J , SMUZ 1 9)  ,SMLZ19)  ,SMGWZ19I  ,FCGK(12)  .CC.N0UZ19) 
0008  DIMENSION    CCN0LZI5)  ,  CONDOM  9)  ,CTP ( 25 ) , ASTAT I  25 ) ,CTR( 25, 4) 

0C39  DIMENSION    AREA (9)  ,  I  ARE  A  (9,2),  I  PLAN!  9)  ,  I  ,\C  PC  P  (  9  )  ,MGS  SP  (9  )  ,MGSEP(9) 

0010  DIMENSION    C.NS  I  17)  , CNS31  7)  ,CNS  51  7)  ,CNS  121  71  ,  FRECI2  5)  ,CTPR(25) 

001  1  DIMENSION    CTRD AY l4),AETUI9),AtTL(91,IAACC(9),IAETI9),NlA19) 

0012  DIMENSION    AMC  MTHI  1  3  )  ,  AVLF  CL  (  1  2)  ,  AVLFCUl  12  )  ,C  (  1 2 )  ,  FCL  (  12  )  ,FCU  (  12  ) 

0013  DIMENSION    KCR0PI7, I  2 ) , NDI Ml  12  ) , PD ACCT (  1 3 ,  11 ) , PREC ( 3  1)  .PSJNSI  12) 
0  014                                      DIMENSION    PWPLZU2)  ,?W°UZ  (12)  ,RA<  12)  ,RCM(  12  ,7)  ,RCN(  12 , 7 ) , RHCI  12) 

0015  DIMENSION     SMACC7 I  1 3 , 3 , 9  I , SM SA TL I  I  2) , SMS AT U I  12)  ,  I A VC 1 3 1 )  ,  TM AX  131 ) 

0016  DIMENSION    T"I HI  31) ,U< 12) .WIND (12) , PR  EC AC ( 25,4 ) ,RUNACCI 2  5, 9) 

0017  DIMENSION    DSRNFF ( 9)  ,  A  I NTE R I  9)  .AAETRS19)  , ACHSCM (9 )  , JSPERC (9 ) 

0  013  DIMENSION    NPNQFI9)  ,NDPERC<9),  S.MSATGI  12)  ,NAMESI  6)  ,  DAY  SCSI  91  ,P'«G(  12) 

0019  CCMPLSX*16     KR0P*l6(7)/"«HtAT'  ."SORGHUM*  ,  •CCHN'  ,'SCY    (JEANS', 
I  'PASTURE ', 'ALFAl  FA' ,  'FALLOW  / 

0020  DATA    CNSl/84. 0,86.0 ,86.0, 86.0 ,80.0,83.0,84.0/ 

0021  DATA    CNS3/81. 0,82. 0,62. 0,62. 0,74. J, 79.0 ,78.0/ 

0022  DATA    CN35 / 7 3 . J , 75 .0  ,  75 . 0,  75  .0 ,6  1 . 0, o9. 0 , 69. C/ 

0023  DATA     CN312/61  .0,65. 0.6*. 0,65. 0,39. 0,55. 0,61  .0/ 

0024  DATA    NAMES/ • I RR. • , • INT. • , «RNFF« , • PERC ' , • AET    ','    SM    ■/ 

0025  DATA    AMONTH/' JAN. • , •FE8.« , • MAR. • , ■ APR.1 , 'WAY    ',  'JUNE ',' JULY • , 
I  'AUG. ',» SEPT' , 'CCT.'  ,  • NOV.'  ,'DEC.'  ,'TOT. '/ 

0026  OATAAS7AT/'>0.0','>0.1,,'>C.2*,,>0.3','>C.4','>0.5,,*>C.6', 

1  '>0.7'  rOO.i'iOO.S"  i'>l.!),i,>l.l,i,>l.2'i,>1.3,i,>l.4,i 

2  •>1.5',,>1.6','>1.7','>1.3','>1.9','>2.J','>3.0',,>4.3«, 

3  '>5.0«  ,'>10.'/ 

002  7                                      DATA  AVLFCU/2.6, I .5,2.5  ,2 .4 ,2 .5 ,2 .6 , 2 .4 , 2  .4 , 2  .  4 , 2  .2 ,  1  .5 ,  1 . 0/ 
0028                                      DATA  AVLFCl /2. 7, 2. 9, 5. 7, 2. 5, 6. 7, 4. 2, 6. 6, 3. 3, 5. 2, 4. I, 4. I, 2. 5/ 

0  02  9  DATA  C/0. 2 ,0.2 ,0.1 77,0. 177,0. 177,0. 177, 0. 159, 0 . 159, C. 138, 0. 138» 

1  0.134,0.131/ 

0030  DATA  FCU/4. 6, 4. 4, 4. 5,4. 6,4.5,4.3,4.0,4.0,3. 8,3.5,2.3,1.7/ 

0031  DATA  FCL/9. 4, 9. 4.,  14. 2, 7.0  ,13.9,9.  I,  13.7,6.8,9.2,7.0,7.0,4.3/ 
0  032  DATA  FCGW/ 7. 0,6. 9, 9. 4, 5. 8 ,9.2  ,6.7  ,3.8,5.4,6.5,5.2  ,4.6,3 .0/ 
0033  DATA  PWPUZ/ 2. 0,2. 9, 2. 0,2. 2, 2.0, 1. 7, 1.6, 1.6, 1.4,1.3,0.8,0.7/ 

3  034  OATA  PhPLZ/6. 7, 6. 5, 8. 5, 4. 5, 7.2,4.9,7.1,3.5,4.0,2.9,2.9, 1.8/ 

0035  DATA    PWG/4. 3, 4. 7, 5. 2, 3. 3,4.6, 3.3,4.3,2. 6,2.7,2.1,1.8, 1.3/ 

0036  OATA    SIMS  AT  U/5.  8, 6. 2,  5.  7,  5. 7, 5. 7,  5. 5,  5. 4,  5. 4,  5.  3,  5. 2,  ■4. 8,  4.  8/ 

0  03  7  DATA    SMSATL/1  1.3,11.5,16.3,8. 2,15.9,10.4,15.4,7.7,13.9, 13.2,13.2, 

I  12.9/ 

0033  DATA    SJSATG/3.8,d.3, U. 0,6. 9, 10. 7, 3. 9, 10. 4, 6. 5, 9. 6, 9. 2, 9. J, 8. 8/ 

0  03  9  DATA    U/ 0.4 7, 0.47,0. 39,0.3  9, C. 39,0.39,0.35,0.3  5  ,0.31 ,0.3  1,0.28, 

1  0.24/ 

3040  DATA    N01K/31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31/ 

0041  OATA    fREO/0.0,0.i,0.2,0.3,0.4,0.5,0.6,C.7,0.a,C.9,1.0,l.l,1.2, 
1  1.3,1.4,1.5.1.0,1.7,1.6,  1.9, 2.0, 3. 0,4. 0,5.0.  10.0/ 

0042  OATA    POT, SMMAXL , SMUZPR , SM PR EV  ,E PR  I M/3 7. 5 , 2. 2 , 3 .25 ,12.60 ,0.5/ 
304-3  CATA    AVAILL,AVAILU,CM,  JSCVOL,  OSCAYS.EVAPLK..H,  IRRSLM, 

I  ?EAK,?«6VDS,SN0*,TPARfcA, 7 PREC, JZSMCH  ,M ASThW, *6T/  16»0.0/ 
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FCRTRAN    IV    G   LEVEL  21  MAIN  DATE    =    7SC45  21/27/40 

0044  OATA    AET,  A ETLZ,AETUZ,B2,0PERC,  EXCESS, I  A,  M  ,M  A  ,,MR  ,PACK,PACKPY  ,  PQVOL  , 
I  PGRC,P0NvrL,Pl,P2,P3,SN0MLT, TDPERC,TRNQF,U ,T2/2 3*0.0/ 

0045  DATA    CT ?0 A Y ,W AT ER , S TR VOL, VOLO  IS , CROP , I3AY ,MM, PLCT , SO  I L/ 4*0. 0 , 5*0/ 

0046  OATA    CTP.CTPR  ,CTRDAY,CTR/25*0. 0,25*0. 0,4*0. 0,100*0,0/ 

0047  NAMEL lST/ALPHA/8KUNTA,3RUNTa,CRCP,GCRM,E,GRCW,0SR.ATE,?AVLU,PCVMAX, 
•RCRCP 

0C4  3  NAMf LIST/BETA/NPLCTS,HMAX,IN0ST,L,LTAREA,MSTAR7,S,STCRM,W,YEN0t 

*YSTART,CARCEC,RCTAYR,CYPASS 

0C49  NAMEL1ST/SEED/M00EL,WPCNT 

C 

q  *****    INPUTS    ***** 

C 

C***  REAO    PROGRAM    IDENTIFIERS 

0050  READ(5,5)       NAME , GF , CITY , AND, S TATE 

0C51  5  FCRMAT(20X,5A4I 

0052  RE\Ct5,  10)    WEATH,  MOO  IF  ,.MUO,  EL 

0053  10  FCRMAT(1X,4A4) 

0054  «E\0(5,15)     IN, FLOW, TC.STGR, AGE, PONO 
0C55                               15  FORMAT  UX,cA4) 

C*«*  READ    LEATHER    MGOIFICATICN,     FEEOLOT,     PCNC    ANO    PLOT    PARAMETERS 

0056  REA0(5,SEE0) 

0C57  RE4C15, ALPHA) 

0058  READ! 5.3ETA) 

0C59  PPEVV9=YSTART 

C***  REAO     THE    MONTHLY    AVERAGE    METEOROLOGICAL    OATA 

0  060  PFACI 5,20)     IP  SUNS! I ) ,RHO(  I)  ,R A I  I )  , WINCI I )  ,MMAT (I)  ,1*1,121 

0061  20  FC*MAT(2X,F2.2,F2.0,F4.2,F3.i,F3.1> 

C***  REAO     INPUTS    FOR    SOIL    TYPE,    CROP    ANO    AREA    FOR    E «CH    PLOT 

CC62  READ(5,25)     ( I  I  AREA! I , J)  , J  =  l ,2 > , 1= 1 ,NPLO I S ) 

0C63  25  FCPMAT(2X,I l,2X, 12) 

0C64  RfcAD(5,30)      ( ARE A ( I)  ,  I =1 ,NPL CT S 3 

0065  30  FCRMATIIX.F6.0) 

C***  RE4D    INPUTS    FOR    PLAN     I MPL EMEN TAT  I CN    ANO    CROP    RCTATICN 

0066  READ(5,35)     I  I P L AN  I  J  I  ,  KOT ATE ( J )  ,  J= I , NP LOT S ) 
0C67                                35  FGRMAT(U,2I1) 

0  0o3  P.EAOI  5,40)     (MGS3PU  )  ,OGSBP<  I)  .MGStPlI  ),  OGSeP{  I  ) ,  1  =  I  ,NPLOTS  ) 

006»  40  FCPMAT(4I3) 

C***  REAO    INPJTS    FCR    CAP.CY'S    ECU  AT  ION    ANO    MUNICIPAL    01  SPCSAL 

0  0  70  IFIDARCEO.Ew.U       REAO (5, 45)  (CCNOUZII)  ,CCNCLZ ( I ) .CCNQGW I  I),HUZ(1), 
IHLZt  r ),HGW(  I  1,  I  =  1,NPLCTS) 

0071  45  FCRMATI3F1C.2.3E10.4) 

0072  IFC5YPASS.GE.3)     READ ( 5, 50 ) ( MUNO I S ( I ) , I* 1 , 12 ) 
0C73                              53  rCP.."AT<Fl0.2) 

C 

C*»*  RUNOFF  CURVE  NUMBERS  FOR  SOIL  AND  CROP  TYPES 

0C74  00  55  K=l, 7 

OC75  P.CNI  1  ,<)=CNS1  IK) 

0076  RCN(2,K)=CNSl(K) 

0077  RCNJ3,K)=CNS3«K1 
0C78  PCN( 4,<)=CNS3( K) 
00T9  PCNI5  ,K)  =  CNS5U) 
OCSO  RCN16  ,K)=CNS5(K1 
0031  SCN( 7,K)=CNS5(K) 

0082  »CN(8,M=CNS5(KJ 

0083  RCN(9,M=CNS5(K) 

0034  RCN(1C,K)=CNS3U) 

0035  RCN(ll,X)=CNS5(K) 

0036  55  RCN( 12,K)=CNS12(K> 
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0037 

0C83 

0039 

0090 

0091 

CC92 

0093 

0094 

OC95 

0C9fc 

0097 

0C93 

00^9 

0  100 

Ol.H 

0102 

0133 

0104 

0105 

0  106 

0107 

010? 

0109 

0110 

0  1 1  I 

0112 

0113 

0114 

0115 

0116 

OUT 

0118 

0119 

0120 

0121 

01?2 

0123 

0124 

0125 

0126 

0127 

0123 

0129 

G   LEVEL      21  MAIN  OATE   -    7SC45  21/27/40 

C***       SIZING    PCND    VOLUME    ROUTINE 
A1=L«W 
A2  =  S*(L«-W) 
A3=4./3.*S«*2 
A4-2  .*A2 
A5=4.*S»*2 
C***      VCLMAX    IS    ThE    MAXIMUM    VOLUME    HELD    3Y    THE    STORAGE    FACILITY 

VJLMAX=l Al»HMAXtA2*HMAX*»2*A3*HMAX**3)/3630. 
C»«.«       P3AREA    IS    THE    OIFECT    RECEIVING    AREA    OF    THE    FACILITY 

PSAREA= (  (W+2.*S*riMAX)*(L»2.«SeHMAX)  )/4  356C. 
C***      CALCULATE    THE    CROP    COEFFICIENTS 
00    60    K=l,rPLOTS 
MG30  =  *GSt!P(K) 
C'SP=0GS3P(K) 
MGSE  =  Mf,SEP  IK) 
JGSE=OGSEP(K) 
CRCP=IAREM  K,  1) 

CALL     CRUPCGUROP,MGSB,DGSB,MGSE,OGSE.KCROP,NOIf  ,MMAT) 
INC&CPtK)  MAREAU.l) 
60    CONTINUE 
C*«*       INITIALIZE    VARIABLES 
00    7C    1=1,25 
00    65    J  =  2,9 

IFtJ.LE.4l     PRECACtl ,J)=0.0 
6  5    RUNACCI  I, J  J =0.0 
70    CONTINUE 

DO    75    J=l , 12 
75    KCR0P(7,J)=0.0 

00    80       II=1,NFLGTS 
T3At-EA=TPAHEA«-AREA(  I  I) 
Till  1=0.0 
ECtll  )  =  0.0 
IAET1 II )=0.0 
IAADDU  H=0.0 
OTSVCLl I! >=0.0 
DAYSOSt  II  )  =  C.O 
OSRf.FFd  I  ]=0.0 
AINTERUI}=0.0 
AAETPSl  I  I  J  =  0.0 
ACHSCM)  II  J  =0.0 
OS3??C! II )=0.0 
SMLZI II )=9.35 
SMUZI II  I =3.25 
SMGV'Zt  II  )=6.30 
30     S«PJ(  II  >=SMLZ(II)+SMUZ(  II  ) 
YEAKS=YENO-YSTART+l 
C 
C«**       PRINT     INPUT     PARAMETERS 

WRITE (6, 85)  NAME, OF, CITY  , AND  .STATE ,YST ART, YENC.WEATH.MOCI F ,MCO, EL 
1  ,  IN,FL'uh,TC,STGR,AGE,PC.NO,LTAREA,  STORM,  L  ,  W ,  S ,  HNAX  ,  VCLMA  X ,  PS  ARE  A 
35  FOR.MATI ' I'  ,10X/////////// 10X, 'STATION:' ,3X,5A4, 10X.I4,'  TO  ',14, 
1////10X,  'MGOEL  :  ,,4A4//18X,6A4//18X,'FEECLGT  AREA  =  •,F6.2,'  ACRE 
2S,////1CX,'SIZE  OF  CRITICAL  EVENT:  • ,F4. 2/ /// 10X , ' PCNO  VARIABLES: 
3'//25x,MA)  CASE  DIMENSION--  *,F7.2,'  FEET  aY«fF7.2t'  FEET'//25X,' 
4(31  SIDE  SLOPE--  RUN:RISE  =  '.F3.0,'  :  l'//25X,MC)  MAXIMUM  UEPT 
5H —  ',F5.2,'  FECT»//25X  ,•  [Dl  ^A<!«JM  PCND  VCLUME —  ',F9.2,'  ACR 
6E-INCHES'//25X,ME)  DIRECT  RECEIVING  AREA  (FOR  PRECIPITATION)  — 
7        >fF8.2t'     ACRES'//// 10X,  'ARE  A    V  AS  I  A3LE  S  :  '  ) 
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0130  00    130    J=1,NPLCTS 

0131  PLAREA=AREA(J) 

0132  CPOP=IAREAC J, I) 

0133  SCIL=lAREA(Jf2) 

0134  RPAVLU=PAVLU 

0135  IF(IPLANU).EO.l)  PAVLU=0.0 

0136  V.RITSI6.90)       J  tPL  AREA,KROP<  CROP)  .  SOIL  ,  DSRAT  E.  P  AVLU 

0137  OO    FCRMAT<//15X,'PL0T     »  ,  I  1//25X,  M  A)     AREA —       ',F6.2,'  ACRESV/25X  ,  • 

1(3)  CROP—  '  ,2A3//25X,  »(C)  SOIL  TYPE—  'tlS.'  (  SCS  SOIL  TYPE)'/ 
2/25X.MD)  IRRIGATION  RATE—  '.F5.2,'  INCHES/CAY  •//25X.MEJ  IRR! 
3CATK5N    MANAGEMENT—  IRRIGATION    3ELGW       ',F5.2,'       FIELD    CAPACITY') 

0133  IFUPLAN(J).cQ.l)    WRITE(6,95) 

0139  IFIIPLANM  J)  .EQ.2)    WRITE<6,100) 

0140  IFIIPLANU  J.ED.3)    WRITE(6,110) 

0141  95    FOR"AT(/25X,' (F)     PLAN    IMPLEMENTED —         RUNCFF') 

0142  100    FCRMAT(/25X,MF)     PLAN    IMPLEMENTED —         RUNCFF    ANU    IRRIGATION') 

0143  HO    FCRMAT(/25X,' <F )     PLAN    IMPLEMENTED —  IRRIGATIC.N') 

0144  IFIROTATEIJJ.EQ.IJ     *R I T E( 6, 12 Oi 

0145  120    FORMAT 1/25X, • (G)    CROP    ROTATION    WITH    —  FALLCW) 

0146  PAVLU  =  RP.\VL'J 
0  147  130    CONTINUE 

C 

C 

C  *»»**    ENTER    YEARLY    LOOP    ***** 

C 

0148  DO  980   flY=l, YEARS 
C 

0149  CO  160  1=1,13 

0150  DC  150  J=l,U 

0151  00  140  K=1,NPLCTS 

0152  140  IFU.LE.3)  SM  ACCT  t  I ,  J  ,. 0  =  0.0 
0133  150  PCACCTI I , J)=0.0 

0  154  160  CONTINUE 

C*»*   ESTABLISH  CSCP  ROTATIONS 

0155  DO    1?C    K=1,NPLCTS 

0156  170     IAREA(K.,U=INCP0P1K) 

0157  IFIR0TAYR.E3.1)    GO    TC    190 

0153  DO    180    J=1,NPL0TS 

0159  IF(IAR6A(  J.D.EQ.l)     GO    TO    130 

0160  IF(RCTATE( JJ.EQ.2)    GO    TO    130 
Olol  IAREA(J,1)=7 

0162  180    CONTINUE 

0163  ROTAYP=l 

0164  GO    TO    200 

0165  190    R0TAYR=2 

C***       INITIALIZE    VARIABLES 

0166  200  IOISOA=0.0 

0167  MAXVOL=0.0 
0  169  LKEVPT=0.0 

0169  VOLIRR=0.0 

0170  VCLCHC    =    CO 

0171  IHNY.GT.l  )    ,MSTART=1 

0172  WRITE(6,210) 

0173  210    FORMAT! «1',46X, •**•**    ANNUAL    SUMMARY    **•**') 
C 

C 

C  «*«««    ENTER    MONTHLY    LOOP    ***** 

C 
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0174  00    840      NM=MS7ART,12 
C 

C*»*  ES7A3LISH    CROP    ROTATIONS    FOR    WHEAT 

0175  00   230    II=1,NPL0TS 

0176  IFIINCROPf II) .NE.l)     GC    TO    230 

0177  IFIROTATECI  n.E-3.2)    GO    TO    230 

0178  IFIR0TAYR.CC.2)    GO    TO    223 

0179  IFINM.GT.MGSEPl  II  J  1     I  AREA ( 1  1 ,  1) »7 
0130  GO    TO    230 

0181  220  IFtNM.LT.HGSBP(II))     I  AREA 1 1  1 , 11  =  7 

0132  IF(NM.GE.MGSBP(II)J     I  AREA ( 1 1, 1)=INCR0P( II ) 

0133  230  CGNTIMJE 

C*»«  READ    ''CNTHLY    METECRCLOG IC  AL    DATA 

0134  240  REACI  1.250, END=1520)    KAN,  ST  INO,  YE  AR  ,MCNTH  ,  (  PREC  <  I  >  ,I=«l,31)  , 

KTMAXIIJ  ,1  =  1,31)  ,  CM  INI  I  )  ,1=1  ,31) 

0135  250  FORMAT!  12,  I  4, 21 2 , 3 IF4. 2 ,62F3. G) 

0136  IFtSTINO.NE.  INCST)    GO    TO    240 
0167  IFt YEAR.LT. YSTART-1900)     GC    TO    240 
0189  IFIYEAR.GT.YFNO-1900)     GO    TO    990 

0139  IFe'CNTH.LT.MSTART.ANO.YEAR.EC.YSTART-^OO)     GO    TO    240 

Cl^O  ACTIPR=0.0 

0191  OSCAY=0 

0192  NDIM<21=28 

019?  IF(r,M.FQ.2.AN0.TMAXC29).LT.900)     N0IM(2)=29 

0194  NOAYS=NDIMJNM) 
C 

c 

C  »»»**    ENTER    CAILY    LOOP    ***** 

C 

0195  CO    800      ND=1,N0AYS 
C 

C***  THE    FCLLCMNG    STATEMENTS    CORRECT       FOR    MISSING    TATA    CN    INPUT    TAPE 

0196  IFITMAX(NO)  .GT. 250.0)    TMAXl NO )=PDT*  100. 0 

0197  IF(T»IN(N0).GT-250.0)     TM I  Nt  NO  )=POT-f  100.0 
0133  IF(PPFCINO)  -GT.99.97)     PREC(ND)  =  0.0 

C*«*  TAVG     IS     THE    AVERAGE    CAILY    AIR     TEMPERATURE,     DEGREE    FAHRENHEIT 

0  19  9  7AVG(N0)=(TMAX(\D)+TMININC)  )/ 2.0- 1 00.0 

C***  SUBROUTINE    HTHRMO    ADJUSTS    ThE    PERC IPI TATICN    RESULTING    FRCM 

C***  SEEDING    CLCUOS 

0200  CALL    V.TRMCC  (PKEC^.D)  ,NM,MCOEL  .WPCNTJ 

C***  THE    FOLLOWING    CARD    EVALUATES    WHETHER    THE    24    HCLR    CESIGN    STORM 

C***  HAS    BEEN    EXCEEDED. 

0201  IF(PRECiriD)  .GE. STORM/1.  14)     WRITEI6.260)     NM.NO , YEAR , PREC (NO) 

0202  260  FCRMATI2CX,  12 , '/•, 12, '/',  12,'       CRITICAL    EVENT    EXCEEDED       ', 

12X.F1C.2,1        INCH    STCPM       '  ) 
C 
C 

C  ***  CALCULATION  Of  POTENTIAL  EVAPQTRANS P IRAT I CN  3Y  "EANS  OF 
C  PENMAN  EQUATION  *** 

C 

0203  R=RCRCP 

C***  THE  FCLLOWING  CARD  CHECKS  FOR  SNOW  COVER 

0204  IF(PACK.GT.O.l)    R=0.70 

C*«*  THE    NEXT    ThG    CARDS    CONVERT    TAVG    TO    ABSOLUTE,    DEGREE    KELVIN 

0  2  05  C E\T=  IT AVG (ND) -32.0)* ICO. C/ 180.0 

0206  A3ST=CENT»273.16 

C***  ES    IS    THE    DAILY    CALCULATED    SATURATED    VAPCR    PRESSURE,     IN    ."ILLI8ARS 

0237  ES=33.9*( I 0.03  738*CENH-C. sO 72 )»*3-0.000  0l5*A3S ( 1. 8*CENT+48> 
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0210 

0211 

0212 

0213 

021* 

0215 

0216 

0217 

0213 

0219 
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1       +0.00136) 

IF(ES.LE.O.C)  ES=0.0 
C**«   ESA  IS  THE  OAlLY  CALCULATED  ACTUAL  VAPOR  PRESSLRE,  IN  MILLIBARS 

E3A---ES*RH0(NM)/100.0 
C***   RN  IS  THE  CALCULATED  OAILY  NET  RAOIATICN,  IN  *!>    CF  WATER 

RN=(l-R)*RAKJM)*(0.22*0.5**PSUNS(NM)  )  -2.  0  l0E-C<;»Ae ST*«*» 

uo.sa*(  i.-?runta-3punto*scrtiesa;  )  )*(o.i*o.9»psuns<nm>) 

IFtF.N.LT.O.O)  RN=0.0 
C***       WlNUU    IS    THE    MC.NTHLY    AVERAGE    WlNDRUN,    MILES/DAY    AT    2    METERS    HEIGHT 

WIN0D=(WIND(NM)*2*)  «0.555 
C***       EA    IS    THE    CC.NVECTIVE    LOSSES,    MM    '«ATER 
CA  =  0.Z6«(E+G.0l»/.INCC)*<ES-ES  A) 
E4LAKE=0.2  6*<  EPRIM+0.G1*WINDD)* IES-ESAJ 
IF(TAVGINO))     270,270,230 
270    DELT4=0.0 
GO    TO    290 
230    OELTA=O.O3<*»TAVG<KD>**0.6  73 
290    GA*MA=1-0ELTA 

C**«       PF7     IS    THE    CALCULATED    DAILY    POTENTIAL    E VAPQ TR ANSP IRAT ICN,    INCHES 
0  220  PET=(  (D£CTA*RN)»<GAMMA*EA)  )/2  5.* 

C 

C***       CSLCULATE    LAKE    AND    BARE    SOIL    EVAPORATION 
0  221  RNSOIL=RN*( ( 1.0-0.2  J) /( 1.0- R)  ) 

0  222  RNLAKE=RN*< (l.0-0.C5)/(  l.C-R)  ) 

0  22  3  PET3S=< <0ELTA*PNSCIL)*(GAMMA*£A))/2  5.* 

022*  LAKEVP=(  (DELTA^fiNLAKE)  M GAMMA *6ALAK£)  )/25.4 

0225  POT=TAVG(ND) 

022?)  IFITAVGIND) .LT. 20.0)     PET=0.0 

0227  IF(TAVGC)O)  .LT. 20.0)     PETBS  =  0.0 

0229  IF(TAVSIND)  .LT.2D.0)    LAKEVP=0.0 

0229  DC  300  MQ=l  ,N»LCTS 

0230  UADO(M3)=IAET(MO)-PET 

0231  CC3GP=[  AREA(MCl) 

0232  IF(KCRCP<CCRGP,NM).EQ.O.O]  IAAOD<  MO  )  =  IA£T  (MQJ-PET  BS 

0233  IF{ lAAOClMC)  .GT.O.l )  IAADDt MO )  =  0.  1 
023*  IF(  lAADD(MC)  .LT.O.O)  I  A  At.  C  (  Mg  )  =  0  .0 

0235  300    CONTINUE 

C***  SUBROUTINE    SNGfcRT    CALCULATES    THE    MOISTURE    AODEO    TO    THE   DISPOSAL 

C»»*  SITE    DUE    TG    SHOWMELT    CN    THE    AREA 

0236  PRECIP=PKECINJ) 

0237  SNOVAP=0.0 
0233  WATE3=PRtC I P 

02:?  CALL  SNOWRT (PRECIP, WATER, PACK , PET, TAVG(NO) .SNOVAP) 

C 

C 

C  **»    EVALUATION    OF    SOIL    MOISTURE    AND    CALCULATICN 

C  GF    ACTUAL    EVAPQTRANSPIRATION    •*• 

C 

02*0  STRVOL=0.0 

02*1  RUNMOS=0.0 

02*2  JJ=0 

02*3  NNN=0 
02**                          330    NNN=NNN*1 

02*5  IF(r;'iN.GT.NPLCTS)    GO    TO    650 

02*6  STRNOF   =0.0 

02*7  JJ=Jjfl 
02*3                            3*0    CROP=  lAREAf JJ, 1 J 

02*9  SaiL=IAREA(JJ,2) 
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OSAREA=AREA(JJ  > 

RAIN= WATER +01 SVCL( JJJ/CSAREA 

IF(OISVOL( JJJ.GT.O.O.AND. PREC IP.LT. 0. 4)     GC    TO    410 

IF(P.AIN.LE.O.O)    GC    TO    400 

CALCULATE    SURFACE    RUNOFF    VOLUME    BY    SCS    METHOD 

IF(KC«C?(CROP,NM!.LE.O.O)     GC    TO    370 

IF{SMUZ(JJ).LT.(PWPUZ(SGIL) ♦0.5*AVLFCU(S01L>)>    GO    TO    350 

IF(SML'ZUJ>.GT.(Pr,PUZ(SOI  L)«-0.3*AVLFCU(SCIL}>  )    GO    TC    360 

GC    TO    330 

••"CJIFY    RUNOFF    CURVE     NUMBER    TO    CONDITION     I    ANTECEDENT    MOISTURE 

RCM( SOIL, CROP) =RCN( SCI L .ChOPl *0.3  9*EXP< 0. 009»RCM SOI L , CROP  J  J 

GC   TO    390 

MGOIFY  RUNOFF  CURVE  NUMBER  TO  CCNDITICN  III  ANTECEDENT  MOISTURE 

P.CM(SG!L,CPCP)=RCNISCIL,CRQP)  *l.9  5«EXP[  -C  .0066  2*RCN  (  SCI  L  .CROP) ) 

GC  TO  390 

IF(  SMUZUJJ.LT. 0.6«FCU(  SOIL))    GO    TO    350 

IFISMUZUJ  ).GT.0.9*FCU(S0IL)I    GO    TO    360 

P.CMI  SCIL.CK0P)=RC.N(S0IL,CP0P) 

SI=10CC.O/PCM( SOIL»CRCP)-10.0 

ER=RAIN-Q.2*SI 

IF(ER.LT.O.C)     GO    TC    410 

RN0F=EP.**2/(RAIN+0.8*SI  ) 

GO    TC    420 

EVALUATE    INTERCEPTICN    STORAGE 

RNCF=C.O 

IA=0.0 

GC    TO    430 

»NQT=0.0 

IA  =  0.  I 

IF(  IA.GT.RAIN)     IA  =  RAIN 

IF((IA*IAA0CUJ)).GE.0.1)     IA=»Q.l-IAADC<JJ> 

EVALUATE    INFILTRATION    INTO    THE    UPPER    ZONE 

Pt*C=PAIN-RNOF-IA 

UZEVAP=0.0 

IF1CARCEQ.EG.2)    GO    TC    435 

SUBPCUTINE    DARCRT    EVALJATES    THE    FLCW    WITHIN   THE    SOIL    PRCFILE 

3Y    APPLYING    ThE    GNE-D  IMENSICNAL    OARCY    ECUATION 

0230  CALL  CA.4CPT  (PER C»FCU(  SCI L ) , SMUZ ( J  J  I ,FCL ( SCI L ) . SMLZl J J ) , FCGW< SOIL ) t 
ISMCWZUJ)  ,OPERC,CCNDUZt  JJ)  ,CONDLZ(JJ)  .CCNCGM  J  J  )  i  HUZ  (J  J  )  tHLZIJJ)  »H 
2G«(JJ  ),PhPUZ(S0l4.)  .PWPLZISOIL  ),PWG( SOIL)  ) 

0231  GO    TO    450 

C***   CALCULATE  PRESENT  STORAGE  AVAlLAbLE  IN  UPPER  ZCNE 

0232  435  SVMAXU  =  0.9*S,'SATU(S0IL)-SMUZt  JJ) 

C**«  EVALUATE  WATE*  CASCADEO  TC  LOWER  ZONE  FCR  STORAGE 

0233  PERCL=PERC-SM«AXU 

0234  IF(PFPX.GT.SMMAXU)     PERCYS VMAXU 

0235  1F(P'?PCL.LT.0.0)    PEP.CL=0.C 

0236  IF    CSMUZC JJ).GT.FCU(SOIL) 1    GO    TO    440 
0267  EXCESS=0.0 

02S3  GO    TC    450 

C***       EVALUATE    GRAVITATIONAL    hATER     IN    UPPER    ZCNE 

0289  440    EXCESS=SMUZ(JJ)-FCU(  SCIL) 
C 

C***  IF  THE  CRCP    IS  OORMANT  CR  ThE  SOIL  LIES  FALLOW,  SOIL 

C***  EVAPORATION  IS  EVALUATED 

0290  '50  IFIKCRCPICPCP.NM)  .LE.O.O)  GO  TO  560 
0  291  T!JJJ=0.0 


0250 

0251 

0252 
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C*** 
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0255 
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0264 
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C*«*       MODIFY    PET    BY    THE    PLANT    CONSUMPTIVE    USE    COEFFICIENT 
AET=KCROP<CRO?,N1)»PET 
IF(°ET.LE.  IAETUJ)  )     AET=0.0 

CHECK    WHETHER    SCIL    MCISTUKE    LIMITS    AET    FRCM   THE  UPPER   ZONE 
IFISMUZIJJ  )-<0.3*(AVLFCUl  SU  IL  )  )  +PWPUZ  t  SC  I  L)  J)     460,  460,4  90 
CALC'JATE    AET    FRCM    THE    'JPPER    ZCNE    WHEN    LIMITED    8Y    SOIL    MOISTURE 
AVAILL=SMUZI  J  J  J-PV.PUM  301  L) 
IF  [AVAIL'J.LF.O.O)    AVAILU=C.O 
AETUZ=0. 7*AET*< AVAILU/(0.3*AVLFCU(S0IL)  )  ) 
EVALUATE    AVAILABLE    WATER     IN    THE    LOWER    ZCNE 
AV/\!LL  =  SMLZt  JJ)-PXPLZ(SCIL) 
IF (AVAl LL.LE.O.O)     AV.\ILL  =  0.0 

CHECK    WHETHER     SOIL    MOISTURE    LIMITS    AET    FRCM    THE    LCWER    ZCNE 
lFISMLZtJJ)-<0.3«<  AVLFCLt  SO  11  )  ) +PWPLZ  IS  C  ID)  )     470, 4 7C  ,  430 
CALCJATE    AFT    FROM    THE     LOWER    ZCNE     -.HEN    LIMITED     eY    SOIL    MOISTURE 
AETLZ=0.3*AET*(AVAlLL/I0. 3* AVLFCL <  SOI D )  ) 
GO    TO    5  00 
AETLZ=AET-AETUZ 
GC    TC    500 

EVALUATE     *.ET    FRCM    9CTH    ZCNES    UNCER    WET    CCNOITICNS 
AETUZ=0.7*AET 
AETLZ=0.3*AET 

AVAIIL=S<LZ( JJJ-PWPLZt  SCI L) 

IFISMLZUJ  >.LE.0.3*(AVLFCLt  SO  IL ) ) +PWPLZ I  SCI D  )    GO    TO    470 
lF(CARCEQ.fcC.l.O)     GC    TO    605 
IFIPEPC-S1MAXU)     510,520,520 

EVALUATE    SCIL    MCISTURE 

SMUZUJ)=SMUZ(JJ)*PERC-A£TUZ-  EXCESS 

SHLZt JJ)=SMLZ( JJ)-AETLZ+tXCESS 

GO    TO    610 

3.'VJZ(JJ)-SMUZ(JJ)*SMMAXU-EXCESS-AETUZ 

SM-<AXl^D.9*FCLtS0lL>-SMLZtJJ) 

IF(?FRCL  +  EXCFSS-SMMAXL)     5 40, 5 40, 5 50 

SMLZl JJ)=SMLZUJ)+PEHCL-A£TLZ+EXCESS 

GO    TO    610 

SMLZUJ)  =  SMLZl  JJI+SMMAXL-AETLZ 

OP^PC^PERCL+CXCESS-SMMAXL 

ZZ    Tn    620 

CALCULATE  EVAPORATION  FROM  BARE  SOIL  SURF  ACE ( SEVAP)  FOR  MONTHS  0CT03ER 

THROUGH  MARCH  CR  kHEN  THE  DISPOSAL  AREA  IS  FALLOW 

AETUZ=0.0 

AETLZ=0.0 

IF(PACK.GT.O.G)  GC  TO  5  90 

IFISMUZIJJ  )  .LT.  tFCUtSOID-UISOID  )  )    GO    TO    570 

EO(JJ)  =  FCUlSOIL)-SVUZUJ) 

!cIS.''UZUj;.CE.FCc<  SOIL))     EC(  JJ)=0.0 

CALCULATE    STAGE    1    SOIL    EVAPORATION 

UZEVAP*PETBS 

EC«  JJ)  =  ECI  JJJ  +-UZEVAP 

lFIEO(Jj).GT.'JISOILJ)     JZ£VAP=EO(JJ)-U(SOIL) 

T( JJ) =0.0 

GO    TC    580 

CALCULATE  STAGE  2  SOIL  EVAPORATION 

T( JJI=T(JJ)+1 

UZEVAP*C(SOIL  >*(T!JJ)**0.5)-C  (SOIL)*!  <T  I  J  J)-l>  **0.5> 

IFIUZEVAP.GT. (PEreS-IAETI JJ ) ) )    UZcVA?=PET SS-IAET I JJ J 

IFIUZEVAP.LT. 0.0)    UZEVAP=C.C 
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0  301 

470 

0302 

0303 

460 
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0324 
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0  33  7  IFISMUZI  JJ)-PWP'JZ(SOIL).LT.UZEVAP)     UZEVAP»SMUZ  < JJ  J-PWPUZ  (SO  ILJ 

0338  GO    TO    600 

0339  590  UZEVAP=0.0 

03*0  600  IF(CARCEQ.EQ.l)    PERC=0.0 

03*1  IF(CARCE<3.EC.l>     EXCESS=0.0 

0  3*2  SMUZt JJ)=SMUZ( JJ)-UZEVAP*PERC -EXCESS 

03*3  IF  (SMUZIJJ  l.LE.PwP'JZlSJIL  >J     SMUZI  JJ  )=PWPUZI  SO  I  L! 

03**  1FICARCEQ.EQ.2)    GC    TC    530 

03*5  605  SMUZI  JJ)=SMUZt  JJJ-AETUZ 

03*6  SMLZUJ)=SML2(  JJI-AETLZ 

03*7  610  IF(S''LZ< JJ) .  LT. PwPLZ < SOIL )  J    A  ETLZ  =  4ET  LZ-<  PWPLZ  (SOIL  I-SMLZU  J  )  ) 

03*8  IFISXLZIJJ  )  .LE.PWPLZI SOIL  )  J     SMLZI J J )  =PWPLZ( SOI L ) 

03*9  IFICARCEQ.EC.l)     GC    TC    620 

0  350  OPFPC=SMLZ< JJ)-G.9*FCL( SOIL) 

0351  IF  ICEPC.LT.O.O)    OPERC=0.0 

0352  IF(SMLZ(JJ}.GT.0.9*FCLIS0IL)I    SML I ( JJ )=0 .9*FCL (SO IL ) 

0353  620  AETUZ=AETUZHJZEVAP 

C*«*  SM  IS  THE  SCIL  MOISTURE  IN  THE  GROWING  ZONE,  IN  INCHES 

035*  Sv( JJ)=SMUZ( JJJ+SMLZ1 JJ) 

0355  [ACT(JJ)»IA*UACDUJJ 

0356  AETU<  jJ)=AETUZ 

0357  AETLf JJI=AETLZ 

0358  NDPERC(JJ)=OPERC 

0359  NIA(JJ)=IA 

0360  NP.-JHF  (  JJ)=RMOF 

0361  IF  <  IPLA.M  JJ)  .LE.2)  STSNOF  =    RNOF 

0362  STRVOL  =  STSVCL  +  STP.NCF*  AREA  ( J  J  J 

0363  GP  TC  330 
036*           650  CONTINUE 

C 

c 

C  ♦**  EVALUATION  CF  VOLUME  USEO  AS  1RRIGATICN  *•* 

C 

C"**  Tl     IS    THE    PREVIOUS    CAY'S    AVERAGE    TEMPERATURE.     IN    FAHRENHEIT 

C«*  OGGREES,    T2    IS    THE    AVERAGE    TEMPERATURE    OF    THE    CAY    TWO    DAYS 

C**»  PPICR    TO    TCCAY 

0365  VCLCIS=0.0 

0366  J0IS3A=0 

0367  NCNT=0 

0363  DO    675    MS  =  I  .NPLCTS 

0369  DISVOLt MS)=0.0 

0370  [F(IPLANlMS).EQ.l)    GO    TO    660 

0371  THi*E0=TAVG(ND)*Tl*T2 
03^2  FRE£ZE»TAVC(ND)*T1 
0373  T2=T1 

037*  T1=TAVG(NDJ 

0375  IF(F?.EEZE.L7.6*.0)     FROZE*  1 

0376  IF(THAWE0.GT.11*.0)     FROZE=0 

C***  WHEN  FROZE  EQUALS  1  THE  SOIL  IS  CCNSIOEREC  TC  2E  FROZEN  IT  IS  THAWEC 

C***  WHEN  FROZE  EQUALS  0 

0377  IF(FROZE.ECl)  GO  TO  660 

C*»*  SUUZ  IS  THE  SOIL  MOISTURE  IN  THE  TOP  12  INCHES  CVER  EACH 

C**«  PLOT;  AVLFCU  IS  THE  AVAILABLE  WATER  CAPACITY  CF  THAT  SOIL. 

C»»*  IRRIGATION  nILL  NCT  CCC'JR  CN  CAYS  THAT  Tht  SCIL  MOISTURE  IS  AT 

C*««  A  LEVEL  GREATER  THAN  THAT  OF  THE  PERCENTAGE  CF  AVAILABLE  WATER 

C*«*  SPECIFIES  ?Y  THE  VARIA3LE  PAVLU. 

0373  S0IL=IAREA(«S,2) 

0379  IF(SMUZIMS)  .GT.(PAVLU*AVLFCU( SCIL) ) ♦? wPUZ i SG I L ) )    GO    TO    660 
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03*0  J01SDA=JDIS0A*1 

0391  NCNT=NCNT+1 

0792  IF(PCNVOL.LT.PCV!MAX*VOLMAX)     GC    TO    660 

0  33  3  CISV0L(MS)=CSRA7E*AREA{,MS) 

C***  IF    THE    PCND    VOLUME    IS    LESS    THAN    THE    VOLUME    REQUIRED   FOR    ONE    FULL 

C***  OAY    OF    IRRIGATION,     IT    WILL    3E    ASSUMED    THAT    NC     IRRIGATION    WILL    OCCUR 

C»«*  ON    THAT    DAY. 

0394  PCNV3L=P0NV0L-DI SVOLt MS) 

0315  IF(PCNVOL. GT.O.O)    GO    TO    6  70 

0  386  DISVnL(MSJ=OISVOL(MS)+PCNVCL 

0387  PGNVCL=0.0 

0338  GC    TC    670 

03d9  660  DISVOL(MSJ*0.0 

0  390  670  SMACCT(NM,3,MS)  =  SMACCTINM,3,MS)*DISV0L(MS)/AREMMS) 

3391  VOL0IS=VCLDIS*0[SV0L<MS) 

0392  IF(OISVOLtMS) .GT.O.O)     DAY SDS( NS)=»OAYSOS I  MS) ♦! .0 

0393  675  CCNT INUE 

0394  IF(NCM.GT.O)     I DI SD A=I D ISC A*JC1 SO A/NCNT 
C***  UPDATF    DISPOSAL    DAY    ACCOUNT 

0395  IFIVCL3IS. GT.O.O)    DSC AY=CS3 AY+1 

0396  IF10YP4SS.NE.1)    GO    TC    680 

0397  RUNOFF=0.0 
0396  GO    TU    730 

C***  SUBROUTINE    STORAG    CALCULATES    ADDITIONAL    LCA0ING    TC    THE    STORAGE 

C***  POisD    DUE    TC    FEEDLCT    RUNOFF    GR    MUNICIPALITY    DISPOSAL 

0399  630  CALL    STORAG ( PI ,P2 ,P3 . PREC I P ,SNOW , FROZE, MCNT H, GRCW , DORM, RUNOFF, MUNO 

1IS(MM),RUNMDS) 
C 
C 

C  **»  CALCULATION  OF  SURFACE  AREA  AND  DETERMI NA I  ICN  OF  SURFACE 
C  EVAPORATION  FROM  STORAGE  FACILITY  **♦ 

C 

C***  THE  FOLLOWING  CALCULATION  EXPRESSES  THE  VOLUME  GF  WATER  IN  THE 

C*»*  STORAGE  FACILITY  FN  CU3IC  FEET. 

0400  730  IF(PCNVOL.LE.O.O)  GO  TO  750 

0401  V=PCNVOL*3A30 

C«**    THE    FOLLOWING    CALCULATIONS    OETERMINE    THE    SURFACE    AREA    OF    Ti-.E    S7CRACE 

C***  FACILITY    AS    A    FUNCTION    CF     STORAGE    VOLUME.       AREA    IS     IN    SCUARE    FEET 

C«*«  VOLUME    IS     IN    CU3IC    FEET.       T He    STORAGE    FACILITY    IS    SHAPED    LIKE    AN    INVERTED 

C***  FRLSTRUM    OF    A    PYRAMID.       INPUT    PARAMETERS    TO    SIZE    THE    FACILITY    ARE    LENGTH 

C***  (L)    OF    THE    BASE    IN    FEET,    WIDTH    CF     THE    bASElW)     IN    FEET    ANO    SLOPE    OF 

C*«*  INSIUE    EMBANKMENTS    GIVEN    AS    A    RATIO    CF    RUN    TU    PISEtS).       IT     IS    ASSUMED 

C**«  THE    POND    DOES    NOT    LFAK.        INPUTS     TO    THE    STORAGE    WILL    <3E    NATURAL 

C*«*  RUNOFF,     FEECLDT    RUNOFF,    MUNICIPALITY    DISPOSAL    ANO    PRECIPITATION. 

C»*«  LOSSES    FROM    THE    PCMO     INCLUDE    EVA^OP.ATICN     ANO    DISPOSAL    VOLUME. 

C»*»  B2    IS    THE    A°EA    OF    THF    SURFACE    LIQUID    IN    SCUARE    FEET. 

0^02  HAPRX=(PGNVCL/VCLMAX)»HMAX 

0403  740  VC=A1*HAPPX+A2*HAPRX»*2*A3*HAPRX**3 

0404  DV=V-VC 

0405  DVDH=Al+A4*HAPRXtA5*HA?RX»*2 

0406  H  =  HA?RX»i)V/UVDH 

C417  IF (ABS(H-HAPRX) .LT.0.1)    GC     iO    750 

0  4C5  HAPFX=H 

0  409  GG    70    7^,0 

0^10  750  IF     (H.CT.hMAX)     H=HM4X 

0*.  1  I  32=  t  W*2.*5*H)*(  L*2.*G«H) 

0412  IF(F3CZE.E0.1)    LAKEVPO.O 

0413  LK.EVPT  =  LK.EVPT»LA,<LVP 


89 


FORTRAN    IV    G   LEVEL       21  MAIN  DATE    *    79045  ■  21/27/40 

041*.  ScVAP=82«(LlKEVP/12) 

C«*«  ScVA°    I-S    THF    VOLUME    CF    WATER    EXTRACTEO    FRCM    THE    STORAGE    FACILITY    8Y 

C»»*  F=EE    SURFACE    EVAPORATION. 

0415  IF( (5EVAP/3630).GT.PCNVCLJ    SE VAP*PCNVCL*3630 

0416  ?C'JVCL  =  PCNVOL-(S5VAP/3630J 

0417  IFtPCWOL.LE.O.O)     PCNV3L=C.O 

r 

C««*  THE    VOLUMES    CF    CALCULATED    RUNCFF,     FEECLCT    RUNOFF,    MUNICIPALITY 

C»»*  DISPOSAL,     AND    PRECIPITATION    FALLING   UN    THE    FACILITY    ARE    ADDED 

C***  TG    THE    VOLUME    CF    WATER     IN    THE     STORAGE    FACILITY     (ACRE-IN). 

0418  ?OHVCL=PCNVCL  HP.U.SDFF«L  TARE  A  J  ♦  (  PR  EC  I  P  *3  SARE  A  J  ♦  STR  VOL*  RUNMDS 
C*»*  THE    VOLUME    CF    WATER    REMAINING    AT    THE    £NC    CF    THE    CAY     IS    EXPRESSED 
C***  IN    ACRE-IN. 

C 

C**«  THE    FOLLOWING    STATEMENTS    DETERMINE    WHETHER    THE    STORAGE    FACILITY    HAS 

Ce*»  OVERFLOWED    AND    IF    SO,     THE    JUANTITY    OISCHARGEC 

0«>19  OSCHRG=0.0 

0420  IF(PCNVOL-VCLMAX)        780,780,760 

0421  760  DSCHRCr=PU\VOL-VDLMAX 
0^22  DSCVCLOSCVCL  +  GSCHRG 

C***  V3LUME    CALCULATIONS    TO    INCREASE    THE    PCNO    SIZE 

0423  CC'JTPL    =    1.0 

0424  PCNTSL    =    CCNTRL*100.0 

0425  VCLCHG    =    CCNTRL*PCNVCL-VOLM AX +VOLCHG 

0426  VOLuXl    =    VOLMAX+VLLCHG 

0427  VC3    =    2.0    *    S    *    HMAX 

0428  VCC  =  < <4./3. >*S»*2)  -  ! VCLMX  1* 3630. /HMAX ) 

0429  VCD  =  VC0**2  -  (4.*VCC) 

0430  VCl  =  SQRT(VCO) 

0431  CI*  =  IVCl-VC8)/2.0 

0432  WRITEI6.770)    M , NO , YF AR , OSCFRG, VCLMXi , PCNTRL » 0  IM 

04^3  770  Fr-.R.^AT</,lX,I2,,/•  ,I2,V'  ,12.  •    -    DISCHARGE   CF    ,,F7.2.'    ACRE-IN   REO 

2UIRES    VOLUME    DF    ',fS.2,'    ACRE-IN    FOR    SF6.2,'     X   CCNTRCL    WHERE    L    » 
?W    =    '  ,F3.2) 

0434  PnNVCL=VCLuAX 

043*  IFIDSCHRG.GE.PEAK)     PEAK  =  DSCHRG 

0436  IF  I  YEAR. GT. PREVYR.CR.CM.LT. I. 0)     MM=MM+l 

0437  ??£Vv^=rCAfi 
0<-33  CM=CM*1.0 

0  4  3°  780  CONTINUE 

C 

C***  UPCATE    SOIL    MCISTURE    ACCOUNT    FOR    EACH    PLCT 

0440  DC    785       KIM.NPLOTS 

0-«4l  SMACCT(NM,2,KH=SMACCT(NM,2.KI  )*PRECIP 

0  442  SMACCT(N«,4,KIJ=SMACCTtNy,4,i<I)«-NIA{KI) 

0  44  3  S  UCCTINM,  5,KI  )=SMACCT(  NM  ,5,K  I  1-t-NRNCF  (KI  ) 

044  4  SMACCTt NM.6  ,K I  )  =  .>MACCT ! NM ,6 ,K I ) *N DP ERC I K  I  ) 

0  44  5  SMACCT(Nrt,7,KI  )  =SMACCT(  NM  ,  7  ,Kl  )  *  AET'JI  KI  )  ♦  AETL I  K I )  -»SNOVA  P 

0  446  SuAC;T!NM,S,K.I  )  =  SMACCTINM,3,M  J+SMIKI  )-S.uFO(KJ  J 

0  44  7  735  SMPC(KX)*SM(KIJ 

C 

C***  UPCATE    PCND    ACCOUNT 

0443  ACriRR=ACTIRP+VCLDI S 

0  449  POACCTINM,  3  )  =  PCACC  T  (  NM  ,  3)  +RUNMOS 

0450  PGACCT(NM,4)^?CACCT(NM,4) 'RUNOFF *LT AREA 

0  451  PDACCT(NM,51*PCACCT(NK,5)  f.TR'vCL 

0452  PCACCTJNM.?.  )  =  >>CACC  T(\M,  8)  ♦SEVAP/36J3 

0t53  PCACCT(NM,9)=PCACCT ( \ * , 9 ) OSCHRG 
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0454  PCACCT(NM,10)=PCACCT(NN,10> ♦< PONVOL-POVCL ) 

0455  IF(MOiEQ.NDAYS)     POACCT t NM  ,  I  1J =H 

0456  POVCL=PONVOL 

045f  IF(PCWCL.GT.MAXVCL)     MAXV CL=PCNVOL 

C 

C***       STATISTICAL    °REC IP  I  TAT ION    ANO    rt'JNOFF    FREQUENCY    CATA 

0456  IF(PRECtND).GT.O.O)     CTP3AY    =    CTPOAYH.O 

0459  IPLOT=NPLCTS 

0460  IF(\;>LCTS.GT.4>    NPL0TS  =  4 

0461  DO    795    11  =  1,25 

0462  IF(P°EUNOJ.GT.FREQI  II)  JCTPII  I )=C T? ( I  I)  +  1  .0 

0463  CNTR=0.0 

0464  00    790    <I  =  1  ,NPLCTS 

0  465  IF(N3NCF(KI  ) . GT .0 .0 .ANO  .  I  I  .  EO  .  1  )     CTRCAYI  KI )  "CTROA  Yt  Kl  J*  1.0 

0466  IFI\P\CF(KI  J.GT.O.O)     CNTR  =  l  .0 

0467  790     IFfNPNOFUI  ).GT.Fu.tQl  II  ))    C  TR  (  I  I  ,  Kl)  =  CTR  {  1 1  ,Kl  )  +  l  .0 

0463  795    IF<PREC(NO).GT.FREQ(II).ANO.CNTR. £3.1.0)    CTPR(  I  U=CTPR( 1 1 }*  1.0 

04o9  NPLCTS-IPLCT 

0470  800    CONTINUE 
C 

C 

C  «*«**    exiT    DAILY    LOOP    ***** 

C 

C***       UPDATE    ACCOUNTS 

0471  PDACCTINM, 1)=AMGNTH(NMJ 

0472  PCACCT(NM,2)  =  <;MACCT  (NM,2t  ll*PSAREA 

0473  PCACCTCK,-1,6>=0SDAY 

0474  P0ACCT(NM,7)=ACTIRR 

0475  CO    810       J=2t10 

0476  810    PCACCTI 13, J)=?CACCT( 13, J) ♦PCACCTI NM, J) 

0477  00    330       ,uP  =  l,NPLCTS 

0478  DC    820       J=2,3 

0479  820    SMACCTI  1  3,  J  ,MP  MSMACC  T(  13  ,  J  ,MP  H-SMACCTl  V*  ,J,MP) 
04  30  SMACCTIMM, 1 ,MP)aAMONTH<NM  ) 

0411  830    S*ACCT(  13,  1  ,:,P)=A,MCNTH(  13  ) 

0432  PCACCT(13,l)=AMONTH(l31 

0433  VCLIRR=VCLIRP*ACTIRP 
043',  840    CONTINUE 

C 

C 

C  **«*«    EXIT    MONTHLY    LOOP    ***** 

C 

0435  OSNC..  =  PACK-PACi<PY 

0486  PACXPY=PACK 

043  7  PC*W=((PC;ACCT(13,2)»PCACCT<13,3)*PCACCT<  13,  4  J +POACCTI  13.  5 J-POACCTl 
113,9) )/(?OACCT( 13,2)+?0ACCT(13,3)*P0ACCT(13,4)+PQACCT(13,5) ))*100. 

043^  WASTWK=rfASTWW+PCwW 

0439  PCACCTI 13,  1 1 ) =PCACCT I  12  ,  11) 

04^0  00    950      KT=1, MULCTS 

0491  OS?NFF(KT)=CSS\FF(KT)  «-SMACCTC  13,5,rtT) 

0  49  2  AINTER<KT)=AINTER(KT)+SMACCTl 13,4  ,KT) 

0493  AAETRS(<T)=AAETRSIKT)*SMACCT(  13, 7,^T} 

0  494  ACHSCIKT  )=ACHSCMCKT)*-SMACCT(  13»8tKTJ 

0495  850    0S?PPC(KT)=0S?ESC(<T)»3"'ACCT{  13.6.K.T) 

04?6  IFRSUN  =  IRRSUM*-VCLlRR 

04V7  TPHEC  =  TPREC«-J,YACCT(13,2,l  ) 

0^99  IF  I  (YEAR+1900)  .EO-YSTART)     ORY  =  SMACC T (  13  ,2  , 1 ) 

0499  IFISMACCTi  13,2,1) .Gc.acT)     « ET =S MACCT (  13, 2 . 1 ) 
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0500  IF(SMACCTU3,2,1).LE.CRY)    ORY=SMACCT [ 13 , 2 , 1 ) 
C 

C***   PRINT  PONO  ACCCUNT 

0501  WRITEt6,360)     YEAR 

0502  860    FDR^ATI »0'  i27X, '  WATER    ACCCUNT    FCR    STORAGE    FACILITY    (IN    ACRE- 

1  INCHES)  -    I*', I2//9X,  ' 

3 '/28X,  'INFLOWS'  ,48X,  ' QUTFL CW S' / 12X , • 

5 '/3X, 'MONTH' ,4X, • PRECI P.  MUNICIPAL         FEECLCT         RUNOFF', 4X,' 

6IRR.     DAYS', 4X, 'DISPOSAL    VOL.  SURFACE    EVAP.  01 SCHARG E' ,4X , • VOL. 

7    CHANGE  HEIGHT' ) 

0503  WRITE (6,370)       ( (POACCT( I , K)  ,K=I  ,  11)  , I »l ,13) 

0  5C4  8  70    FURMATt4X,A4,FU.l,Fl2.1,  if  10  .  I  ,F  10.  J  ,F  It  .  I  ,F  1  7.  1  ,F  13.1  ,F15.  I  ,F  10. 

12) 
C 
C*«*       PRINT    SOIL    MOISTURE    ACCOUNTS 

0505  00    910       JM=1,MPLCTS 

0506  CPCP=  [AREAI J",  I) 

0507  S0IL=IAKEA(JM,2) 

0508  "RITE (6,o30)        JM , KRCP I CRC P >  ,S 01 L , ARE  A < J M ) 

0509  830    F0RMAT(////,60X, 'PLOT    NO.  • ,  13  ,/// /  ,25X, ' CROP — •  ,2A8  ,5X, • SCI L    TYPE- 

!-'  ,13, 5X,  •  DISPOSAL    AREA—  •  ,  F6  .2  ,  •     ACRES') 

0510  KRITE(6,89C)        YE\R 

0511  890  FORMAT  CO' ,35X,  "WATER  BALANCE  (INCHES)  IN  THE  DISPOSAL  AREA  -  19', 

2 V32X, 

3     •  INPUTS'  ,3SX, 'CUTPUTS'/21X,' ',3X,  • 

i, V9X,  'MONTH' , 

57 X, 'PREC IP  I  TAT ICN',4X,  '  IRRIGATION'  ,3X,' INTERCEPT  I CN' ,2X.' SURF  ACE    R 
6UNCFF'  ,3X, 'PERCOLATION' ,3X,  'AET',3X, 'CHANGE     IN     SM • ) 
*RITE(6,900)       ( (SMACCT1 I ,K,JM> ,<= 1,3)  ,1=1  ,13) 
900    FC?MAT( 10X.A4.7F15.2) 
910    CONTINUE 

WR:TE(6,920)        PCViW 
920    FCRMATl'O' ,10X,' PERCENT    CF    WASTEWATER    CONTROLLED* • ,F 1C. 2) 

WRITE     16,930)        IOISCA 
930    FCPMAT  I  'V  ,  10X,  'POTENTIAL    DISPOSAL    DAYS  =  *  ,1*1 

rSItE(6,S4G)        PACK.OSNOW 
940    FCRMATI  '0'  ,  10X,  'PACK.    CN    0ECEM8ER    31    =',F5.2,15X, 
I 'CHANGE    IN    SNCrt    STCRAGE=' ,F5.  2) 
WRITE(6,950) 
950    FCRMATI  '0'  ,  10X,  •  INPUT  S-CU  TPUT  S-CHANGE    IN    SNCW    STCRAGE'CHANGE    IN 
1 S  T I L    MOISTURE'  ) 
M\XV0L  =  '',AXVGL*1CC.0/VCLMAX 
w3IT=(6,36C)       MAXVCL 
960    FCRMATI »0« , 10X, 'PERCENT    GF    MAXIMUM    POND    VOLUME    REQUIRED    -*»F7.2J 
EVAPLK=EVAPLK+LKEVPT 
hRITE(6,97C)       LKEVPT 
970    FfS"AT(  '0'  ,  10X, 'ESTIMATED    LAKE    EV  APOP.AT  ICN ,     INCHES    =',F6.2) 
980    CONTINUE 
C 
C 

C  *****    EXIT    YEARLY    LOOP    *•*** 

C 

0530  990  CONTINUE 

0531  1520  CONTINUE 

C***   CALCULATE  AVERAGE  ANNUAL  VALUES 
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0532  EVAP=EVAPLK/YEARS 

0533  CMNEW=CM 

0534  IFIMM.EO.OJ       MM»l 

0535  CCUNT=CM/MM 

0526  IF(CCLNT.EQ.O.C)    M«=0 

0537  IF(CM.EC.O.O)    CM^YEARS 

0533  OSCRG=OSCVCL/CM 

0  539  CM=CMNEH 

0540  CCNTP.L  =  WASTWW/YEARS 

0541  IP.avCL=IRRSUM/YEARS 
0  542  APREC=TPREC/YtARS 
05't3  RANGE*WET-DRY 

0544  AVGMO=FVA?-APREC 

0545  00    lCOO    J=1,NPLCTS 

0546  DS?EKCU)=OSPERC(J)/Y£ARS 

0547  CSRNFFU)  =  CSRNFF(J)/YEARS 

0548  ACHSCft J)=ACHSCV( Jl/YEARS 
0  54)  AI'iTERl  J)  =  AINTER(  J)  /YEARS 
0550  AAETRSt J ) = A AET kS (J 1/YEARS 
0:51  1000  OAYSOS(JJ=GAYSDS(J)/YE>RS 
05c2  IFl JPLuTS.CT.4)  MPLCTS=4 

0553  DO  1020  J=l  ,25 

0554  PRECACt J,3)=CTP(J> 

0555  IF(CTPGAY.CT.G.O)    C TP ( J )=CTP( J)/CTPOAY* IOC. 0 

0556  PRECACt J,l)=ASTAT(J) 

0557  PRECAC(J,2)=C7P( J) 

0553  PRECACI J,4)=CTFR(J) 
05c9                                        DC    1010    1  =  1  ,i\PLCTS 

0560  INUMI=5+I 

0561  RUNACCI J,IMJM11-CTR( J,  I) 

0562  IFICTROAYI  I  J.GT.C.O)    CTRI  J ,  I )  *CTR  (  J  ,  I  )/  CTRO  AY  (  I}*  100.0 
0  56  3  INLW«I*1 

0554  1010    RUNACC(JtINLIM)=CTRUt!l 

0565  1020    RUMACC! J, 1)=ASTAT( J) 

0566  NPLCTS=1PLCT 
C 

C*«*       PRINT     INPUT    PARAMETERS 
0  567  wRlTE<6,35  )       NAME, OF , CITY , AND .STATE, YST ART , YENO ,WcATH  ,MCDIf . MCOrSL 

l.IN.FLCW.TCSTCR.AGE.FCNO.LTAREA.STORhSL.h.S.hlMX.VCLMAX.PSAREA 

0558  30    10H0    J=1,NPLCTS 

0569  CROP=iNCi<CP(J) 

0570  PLAP.f  A=APEA(J) 

0571  S0IL=IAREA(J,2) 

0572  R?AVLU  =  PAVL'J 

OS-*!  IF(IPLAN(J).E0.11    PAVLU=0.0 

0  5  74  WRITE (6 , 10 30 ) J , PLAR E A, KROP ( CROP ! , 30  IL  .DSRATE ,PAVLU 

0575  1030    Fnp."\T(//15X,'PLCT     •  ,  U//  25X,  •  ( A  J     AREA —       ',F6.2,'         ACRES  '//25X,  • 

1(d)  CRCP--  ',2A8//25X,*(C)  SOIL  TYPE—  ',13,'  (SCS  SOIL  TYPE)'/ 
2/25x,'(D)  IRRIGATION  RATE—  «,F5.2,'  INCHES/CAY  »//25X,ME)  IRRI 
3GATIGN    MANAGEMENT—  IRRIGATION    3EL0H       '  ,F5.2,*       FIELl)    CAPACITY') 

IFUPlANIJl.EQ.l)  WRITE! 6,1040 
IF(  I  PLAN! J) .ci.2)  WRITE(6,105  0> 
I c  t  i  ? L AN ( J ) . F Z . 3 )     «RI TE  < 6 , t C6C) 

F0RMATl/25Xf  •  (F)    PLAN    I  IMPLEMENTED—         RUNOFF') 

FORMATI/25X,' (F)     PLAN    IMPLEMENTED —         RUNCFF    Af*D    IRRIGATION') 
FCR*ATl/25X,' IF)     PLAN     IMPLEMENTED—  I  RP  IGAT  1  C.\«  ) 

IF (ROTATE (J). EC. I]    WR IT E I  6,  10 70 ) 
F0PMATI/25X,' (G)    CROP    ROTATION    WITH    —         FUl.Ck'1 


0576 

0577 

0578 

0579 

10*0 

0530 

1050 

0  531 

1060 

0  5  32 

0533 

10  70 
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05*4 

0  53  5 

0536 

0537 

0538 

0589 

0590 

05<n 

05^2 

0593 

0594 

0595 

0596 

0597 

0598 

0599 

0  600 

0  601 

0  602 

0603 

06C4 

0605 

0606 

0637 

0603 

0609 

0610 

0611 

0612 

0  613 

0614 

0615 

0616 

0617 

0618 

0619 

0620 

0621 

0622 

0  6?  3 

0624 

0625 

0626 

0627 

062? 

0629 

C630 

0631 

C632 

PAVLU=RPAVLU 
1080    CONTINUE 
C 

C**«       PRINT    FINAL    SUMMARY 
WRITE (6, 1090) 
1090    FORMAT!////, 47X, ******    FINAL    SUMMARY    *»**♦•) 
WRITE (6,1 100) 


1100    FORMAT t 

f>KlTE(6 
1110    FORMAT! 

WRITE  (6 
1120  FCPMAT! 

WRITE  (6 

1130  FORMAT ( 

10F',F6. 

WRITC  (6 
1  140  FOf-MATI 

WRITE (6 
1150  FORMAT! 

WRITF<6 
1160  FORMAT! 

WRITE (6 
1170  FORMAT! 


IE  =  ' ,F6.2) 


WRITEI6 
1180  FORMAT! 

WRITEI6 
1190  FORMAT! 

WRITE  (6 
1200  FCSWATt 

WRITE  (6 
1210  FCRMATI 

wP.lTc<6 
1220  FORMAT! 

WRITEI6 

1230  FORMAT! 

li«  ACRE 

WRITEI6 
1240  FORMAT! 

00  1320 

WRITEI6 
1250  FURMATl 

WRITE  (6 
1260  FORMAT! 
IS*  ) 

WRITE (6 

1270  FORMAT! 

1  INCHES' 

WRIT; (6 
1230  FORMAT! 

WRITE  16 

1290  FORMAT! 

I  INCHES 

WRITEIC 
1300  FORMAT! 
1F6.2,' 

WRITE (6 
1310  FORMAT! 


0'  ,10X, 'METEOROLOGICAL  SUMMARY'  ) 

1110)  EVAP 

0',25X, 'AVERAGE  ANNUAL  LAKE  E VAPORATl CN*' ,F6.2, '  INCHES') 

1120)  APREC 

0' ,25X, 'AVERAGE  ANNLAL  PRECI P I T AT ICN='  ,F6 .2 , '  INCHES') 

11 30  >  PANGS, DRY, WET 

0'  ,25X,'?RECIP!TATICN  RANGE' • , F6 .2. '  INCHES    (FROM  A  LOW 

,'  INCHES   TO  A  HIGH  OF  ',F6.2,'  INCHES)') 

1140)  AVGMO 

0'  ,25X, 'AVERAGE  ANNUAL  MOISTURE  OEF  IC  I  T* »  ,  F6.2, '  INCHES') 

1150) 

0'  ,  10X,' SUMMARY  OF  FONO  OPERATIONS') 

1160)  MM 

0',25X,'NO.  OF  YEARS  HAVING  A  0  I  SCHARGE*' ,1 6) 

117C)  COUNT 

0' ,25X, 'AVERAGE  ,NC.  OF  DISCHARGES  /  YE*R  HAVING  A  OISCHARG 


118C)  OSCRG 

0'  ,25X, 'AVERAGE    D ISCHARGE= ' ,F6. 2 ,  IX , ' ACRE- 1 NCHE S' ) 

1190)     CCNTRL 

0' ,25X, 'AVERAGE  PERCENT  OF  WASTEWATER  CONTROLLED"' ,F6.2) 

1200)  CSCVOL 

0'  ,25X, 'TOTAL    OISCHARGE    VOLUME* ', F9 .2 , •    ACRE-  INCHES • ) 

1210)     CM 

0'  ,25X. 'TOTAL    NO.    OF    DISCHARGES* • iF4.0l 

1220)     PEAK 

0'  ,25X, 'MAXIMUM  0  ,' SCH  ARGE=  »  ,F  6.  2  ,  '  ACRE-INCHES') 

1230)  IhRVOL 

0* ,25X, 'AVFRAGE  ANNLAL  VOLUME  OF  WASTEWATER  APPLIED*' ,F8.2 

INCHES' ) 

1240) 

0' ,  10X, 'SUMMARY  OF  DISPOSAL  PLOTS') 

J=l,;-iPLCTS 

1250)  J 

C  ,  15X,  'PLOT  ',  II  ) 

1260)     DSRNFFU) 

0' , 25X, 'AVERAGE    ANNUAL    DISPOSAL    AREA    RUNOFF*' ,F6.2 , •     INCHE 

127:)  OS?FRC(J) 

0' ,25X, 'AVERAGE  ANNUAL  DISPOSAL  AREA  P ERCCLATIQN-* ,F6 .2 , ' 

1230)  DAYSOS(J) 

0'  ,25X,' AVERAGE  ANNUAL  NO.  OF  DISPOSAL  0AYS-= '  ,F6.  1) 

1290)  AINTERJJI 

0» ,25X, 'AVERAGE    ANNUAL    DISPOSAL    AREA    I  NTS RCEPTI CN  =  '  ,F6.2 , • 

) 

130C)     AAETRSIJ) 

0' ,25X, 'AVERAGE    ANNUAL    DISPOSAL    AREA    E VAPCTRANSP (RATION*' , 

NCHES* ) 

1310)    ACHSOM(J) 

0'  ,25X, ' AVERAGE  ANNUAL  DISPOSAL  AREA  CHANGE  IN  SOIL  MCISTU 
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IRE=',F6.2,'     INCHES') 

0633  1320  CCNTINUE 

0634  WRITE  16,1330) 

0635  1330  FORMAT! '0« ,10X, 'SUMMARY  CF  STATISTICAL  CATA  '  ) 

0636  WRIT=(6, 13A0) 

0637  1340    F^RMATCO' ,41X, 'PRECIPITATION    FRE3UENCY    CATA ', //27X,  •  INTENSITY' ,  5X 

I, 'FREQUENCY',  5X,« FREQUENCY'  ,5X,'RUNCFF    FRfcQ  .  '  ,/29X,  •  i  IN.  )  • ,  I  OX,  •  U 

2),,9X,'{0AY3)',10X,'(CAYS)',/) 
0  63  8  V.  RITE  (6, 1350)  ({  PR  EC  AC  t  I, J  ),J=  1, 4),  1  =  1, 25) 

0639  1350    FGRMAT(2C'X,A4,3F15.2) 

06^-0  WRITE(6,1360) 

0641  1360    FCPMATl '0'  ,///60X, 'RUNOFF     FREQUENCY    C AT  A • ,//2 7X, '  INTENS I TY • , 15X. ' F 

IREQUFMCY    (  %  )  •  ,26X  ,  •  FREQ'JE  NC  Y     (C AY S> • , /29X  ,'  < IN  .  >•  ,7X, • PLCT    1  PLO 

3T    2         PLOT     3         PLCT    4',3X,'PLCT    I  PLOT    2         PLOT    3         PLOT    4',/) 

0  6<-2  WRITE  (6, 13  70)  (  (RUNACCt  I,J),J=l,9),I=l,25) 

06'»3  1370    FCRMAT129X.A4,  5X,4F9.2,4X,4F9-2) 

0644  STCP 

0645  END 
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0001  SUBROUTINE  CKCPCO  (  CPCP  ,MCSB,  CGS8  ,MGS  fc  ,  CCSE  ,KCPCP  ,N0[  N,  MMAT  ) 
"***  SUBROUTINE  CPCPC3  CALCULATES  THE  CROP  COEFFICIENTS  FOR  USE  IN 
C*«*  THE  MAIN  PROGRAM.   THE  CR.CP  COEFFICIENTS  ARE  CALCULATED  BY  THE 
C*«*  PROCEDURES  OUTLINES  IN  THECHNICAL  RELEASE  UG    21,  IRRIGATION 
C***  WATrP  REQUIREMENTS,  UNITED  STATES  DEPARTMENT  OF  AGRICULTURE, 
C***  SOIL  CONSERVATION  SERVICE,  ENGINEERING  DIVISION,  APRIL  1S67. 
C***  SLIGHT  MODIFICATIONS  HAVE  3EEN  MADE  FCR  ACAPTAICN  TC  THE  1CDEL. 
C***  EQUATIONS  FOR  THE  CROP  GROWTH  STAGE  COEFFICIENT  CURVES  *ERE 
C***  DEVELOPED  WHICH  ELIMINATES  THE  NECESSITY  CF  REACING  ThE  VALUES 
C***  FR'j"  THE  CURVES.   INPUTS  TO  THE  SUBROUTINE  INCLUDE  THE  CROP, 
C***  MONTH  AND  CAY  GROWING  BEGINS  AND  ENOS,  NUMBER  CF  CAYS  IN  EACH 
C»*»  MONTH,  AND  THE  MEAN  MONTHLY  AVERAGE  TEMPERATURES  IN  FAHRENHEIT 
C***  DEGREES. 

0002  INTEGER  C RCP  ,  OGSB ,DGSE 

0003  INTFGER  NO  I M(  L2 ) , SH I F T 

0004  REAL  MIDI  12 > / 1 2*0. / , CBMDt  12 ) /  12*0 ./ , ACC (  1 2 ) / 12 *Q. / . PCGS I 12)/ 12*0./ 

0005  REAL  ,'|MAT(  12)  ,«T{12)  ,KCRCP(  7,  12  )  ,PCGS1(  12  ) 

C**«  MGSD=    HONTH    GROWING    SEASCN    BEGINS    EXPRESSED    NUf-ERICALLY    IS    1-12 

C***  DGS3  =    DAY    GROWING    SEASCN    BEGINS    EXPRESSEC    NUMERICALLY 

C***  MGSE=    MONTH    GROWING    SEASON    ENDS    EXPRESSEC    NUMERICALLY    IE    1-12 

C*»*  DGSE=    CAY    GROWING    SEASON    ENCS    EXPRESSED    NUMERICALLY 

C**«  Mio=MEDIAN    GATES    CF    THE    MCN1HS    IN    THE    GROWING    SEASON 

C***  DBMO=    CAYS    BETWEEN    MID    CATES 

C***  4CC=    ACCUMULATIVE    CAYS     IN    GROWING    SEASON 

C***  PCGS=    PERCENT    OF    GROWING    SEASCN    REAChED    AT    HID    CATES 

C*«*  MMAT=MEAN    MONTHLY    AVERAGE    TEMPERATURES 

C*»»  M0S?I=TE-1P0RARY    STORAGE    FCR    MGS8 

C***  MGSF1=TEMPCRARY    STORAGE    FOR    MGSE 

C***  PCGSi=TEMPORARY    STORAGE    FOR    PCGS 

0006  H&SEl=MGSB 

0007  M3SE1=MGSE 

0008  IFtMGS3.GT.MGSE)    GO    TO    10 

0009  GO    TC    20 

C***  WHEN    MGS3     !S    GREATER    THAN    MGSE    SUCH    AS    IN    WINTER    WHEAT    THE 

C***  SUBROUTINE     "SHIFTS"    CR    AJCS    1    TC    MCSB    ANC    MGSE    UNTIL    PGS3    =     13 

C*»*  WHICH    CORRESPONDS    TO    JANUARY.        THIS    SHIFT    WAS    NECESSARY    TO 

C***  FACILITATE     PRCGRAW    LOCPING.       AFTER    CALCULATIONS    ARE    MADE    THE 

C**«  CROP    COEFFICIENTS    ARE       "SHIFTED"    BACK    TC    THEIR    ORIGINAL    MONTHS. 

£**««*»*^«*«*«**»******«********«**«***«***  **a****»*«**********«**«*a»e**3«**« 
^txitottinnj.nntdtJutMMMmuttuuntt.  *«si»****»*»********a«***««i>*v** 
C***  *** 

c**»  **» 

c*»*  *** 

C**«  ***    CAUTION    TO    USER     ***       THIS    ROUTINE    WILL    NOT    WORK    IF    ThE  **» 

C***  GROWING    SEASON    EXCEEDS    ONE    YEAR.  *** 

C**«  *»* 

C***  *** 

c««*  »»» 

,'-_43ttstJ»».t«*M,*»e»*«»«i>?«»*»«**<t«*»»»«»*»»>*«j«i,««n»,»««»jni»«»*f.»i;'cj.>i« 

0010  10  SHIFT=13-MGSB 
OCU  MGSE  =  *GSE«-SH!FT 
0012  MGS3=1 

0C13  20  NPLUS=MGS8+1 

0014  NMINUS=MGS£-1 

0015  MID(MGS8)»< (NDIM(MGSa)-DGS3)/2.)+0GS3 
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0016  DC  30  N=NPLUS,NMINUS 

0017  30  MI0(N)=NDIM(N»/2.O 
0013  MI0(MGSE)=DGSE/2.0 

001<5  03M0(  wGS3)  =  >'I0(MJSB)-0GSa 

0020  CC    40    N=NPLUS,NMINUS 

0021  40    DBM0CN1 =NDI M( N-l ) -M I D < N-l )»PIC(N) 

0  02  2  CBMn(MGSE)=NDIM(MGSE-l)-MlOlMGSE-l)*OGSE 

0023  ACCtPGSB) =C3MC(MGSB) 

00?4  00    50    N=NPLLS,MGSE 

0025  50    ACC(NJ  =  ACC(N-l  J*0tlMDtN) 

0  026  ACCMGSE)=ACC(yGSE)-uI3(MGS£) 

0027  DC    60    ?:  =  HGS3,«GSE 

0C28  60    PCGSINJ*(ACC«N)*100.)/(ACC(MGSE)*MIO(MGSE1J 

0029  IFIMGS81.LE.MGSEU    GO    TO    100 

0030  CO   80    N=l,  12 

0031  NN*N-SHIFT 

0032  IFINN.LE.01    f.N  =  NN*12 

0033  IFINN.GT.MGSEl.ANO.NN.LT.HGSBH    GO    TO    70 

0034  PCf,Sl(NN)=PCGS(N) 

do:?  gc  tc   so 

0036  70    °CCSIU,N)=0.0 

0037  30    CONTINUE 

3  33<<  DO   90   N= 1,12 

0C39  90    PCGSINI=?CGS1 (N> 

0040  100    MGSB=*GS3  1 

0041  MCSE=MGSEl 

0042  DO  110  J  =  l  ,12 

C***   <T  IS  A  CLIMATIC  COEFFICIENT  APPLIED  TO  THE  CRCP  GRCHTH 
C*«*   COEFFICIENT.   IT  IS  CALCULATED  3Y  THE  FOLLOWING  ECUATICN: 

3C43  KT(  J)=.0173*,«"AT(  J)-. 314 

0044  IFIMMATU).LT.36.)    KT(J>».3 

0045  113    CCNTIMJE 

C***       CR0?=1    FCR     WHEAT 
C**e       CRCP=2    FOR    SORGHUM 
C***      CPCP=3    FCR    CCRN 
C***      CR0?=4    FCR     SCY?EANS 
C**«       CA0?=5    FCR    PASTURE 
C«*«       C:>Z?  =  b    FOR    ALFiLFA 
C***      CP3r,=  7    FOR    FALLCW 

0046  GC    TL    ( 120,130,  140,150,  160,  130,200) .CROP 

0047  120    X0AR=50. 

004S  A= 1.39093399                                                                                                                                                                  WHEAT 

0049  3=-0. 00368373  WHEAT 

0050  C=-0. 00004976  WHEAT 

0051  D=-0. 00000233  WHEAT 

0052  c=-C.C3000C04  WHEAT 

0053  GC  TC  213 
0C54  130  XeAR=50. 

0055  A= L.0552i355  SORGHUM 

0056  3O.CC19d6C0  SURGHJM 

0057  C=-0.C0051577  SORGHUM 

0058  D=C.000P0C45  SORGHUM 
00*9  5=3.00.100011                                                         SORGHUM 

0050  cc  rn  210 

C061  140    Xa4P=50. 

0C*>2  A=L.  02305323                                                                                                                                                           CCRN 

0C63  B=0.0C880046                                                                                                                                                           CORN 

0064  C=-0. 00031919                                                                                                                                                        CORN 
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0065  0=-0.C0000194                                                                                                                                                       CORN 

0C66  E=). 00000007                                                                                                                                                          CORN 

0067  GO    TO    2  10 

0068  150    XSAR=50. 

0C69  A=0.7479C430                                                                                                                                               SOYBEANS 

0070  3=0.01474796                                                                                                                                                         SCY8EANS 

0C71  C  =  -0. 00013486                                                                                                                                                      SOYBEANS 

0072  0=-O.CC0GC443  SOYBEANS 

0073  E  =  0.  SCY3EANS 

0074  CO  TO  210 

C***       POP.    PERENNIAL    CROPS    SUCH    AS    ALFALFA    ANO    PASTURE,    VALUES    OF    THE 
C***       CROP    CCEFFICIENTS    ARE    3EST    PLOTTED    ON    A    MONTHLY    BASIS    THEREFORE 
C***      EQUATIONS    k^ERE    NOT    OEVELOPEO.       MONTHLY    VALUES    HERE    INTEGRATED 
C***       WITHIN    THE    ROUTINE    FOR    PASTLRE    ANO    ALFALFA. 

0C75  160    KCRCP(5,1 1=0.49 

0076  KCR0P(5,2)=C57  PASTURE 

0077  KCRCP15,3)=0.73  PASTURE 
0073  KCRPP{5,4)=0.85  PASTURE 
0079  KCRnP(5,5)=0.90  PASTURE 
3030  KC^CP(5,6)=0.92  PASTURE 
0081  KCRCPI5, 7)=0.92  PASTURE 
3C32  KCRUP  (5,81=0.91  PASTURE 
0CS3  KCRCP15,9)=C.o7  PASTURE 
0034  KCRC^(5, 10)=0.79  PASTURE 
0C85  KCRCPI5, 11 1=0.67  PASTURE 
C036  KCRCP(5,12>=0.55  PASTURE 
0037  OQ    1/0    J=L,  12 

0033  KCRt^?(5  ,J)=KCRCP(5,J>  *KT(  J)                                                                                                            PASTURE 

0C39  IF<PCGS( Jl.LE.O.O)     KCROP(  5,  J)  =0. 0                                                                                           PASTURE 

0090  170  CONTINUE 

0091  GO  tc  230 

0092  130    KCPCP16, 11=0.63 

0C93  KC3CP<6,2)=0.73                                                                                                                                   ALFALFA 

0094  KCRfP(6,31 =C.  S6  ALFALFA 

0095  KCRGP  16,41  =  0.99  ALFALFA 
J0J6  K.CR0P(6,5  1  =  1.  J8  ALFALFA 
0097  KCRC°(6,6)=1.13  ALFALFA 
0C93  KCR0P16, 71=1.11  ALFALFA 
0C99  KCRPP(6,8) =1.06                                                                                                                                                   ALFALFA 

0100  KCRCP(6,9)=C99  ALFALFA 

0101  KCRCPI6, 101=0.91  ALFALFA 

0102  KCRCPJ6, 111=0. 78  ALFALFA 

0103  KCKGP16, 12>=0.64  ALFALFA 

0104  DC    190    J  =  1  ,12 

01C5  <CR0P(6,J)=KCRGP<6,J)*KT( J)                                                                                                                  ALFALFA 

0106  190    IF<PCCS< J)  .LE.0.01    KCROP( 6, J) »0.0 

01C7  GO    TC    230 

0108  200    X3AP=0. 

0109  A=0.  FALLOT 

0110  3=0.  FALLOW 

0111  C=0.  FALLOh 

0112  0=0.  FALLOW 

0113  E  =  0.  FALLOW 

0114  210    00   220    J=l,12 

0115  Z=PCGSl J1-X3AR 

0  116  KCROP(CROP,J)= (A*B*:*C*Z**2+0*Z**3*E*Z*»4)«KTt J) 

0117  IF  (PCGSI Jl.LE.O.O)  KCRCP (CROP, J) =0 .0 

0  118  2?0  CONTINUE 
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FORTRAN 

IV 

G 

LEVEL 

0119 

230 

c»»* 
c*»* 

0120 

0121 

0122 

0123 

21  CROPCC  OATE  =  7S0A5  21/27/40 

CCNTINUE 

SINCE    THE    MAIN    PROGRAM    APPLIES    THE    CRCP    COEFFICIENT    (KCROP)     TO 

THE    ENTIRE    .MONTH,    THE    KCRCP    WAS    PROPORTIONED    ACCORDINGLY    TO 

COMPFNSATE     FOR    THIS.        THE    NEXT    TWO    CAROS    CO    THIS. 

KCROP(CR0P,MGSe)=KCRCPt  CROP  .MGS3)  *<  NO  IM  <  ,«GSa  1-OGSB*  1)  /N0IMJMGS31 

<CRCP{CRCP,yGSE)=KCRCP(CRCP,MCSE)*OGSE/NDIM<MGSE) 

RETURN 

END 
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FORTRAN    IV    G    LEVEL  21  WTRMQO  DATE    =    7SC45  21/27/40 

OOOl  SUBROUTINE    WTRMCC ( PREC I P, f ONT H, MOCEL, WPCNT) 

C**»  SUBROUTINE    WTRMOO    ADJUSTS    THE    PRECIPITATION    RESULTING    FRCM 

C»»*  SEEDING    CLOUDS 
C 

C«»*  NO    PRECIPITATION    MOOIFICAIICN 

C*»*  .MODEL    2:       INCREASE    PRECIPITATION    DURING    MARCH    TFRCUGH    SEPTEH3ER    3Y 
C***  THE    FOLLOWING    PERCENTAGES    -    755    FOR    RAINFALLS    <    3-10    I  .M , 

C*»«  30?    FOR    RAINFALLS    <    0.50     IN,     10?    FOR    RAINFALLS    <     1.0    IN, 

C**«  -10"?     FOR    RAINFALLS    >     1.0     IN 

C*«*  MCOEL    3:        INCPEASE    ALL    RAINFALLS    DURING    EVERY    fCN  TH    3Y    A    SPECIFIED 
C*«»  PERCENTAGE,     rfPCENT 

C*»*  MODEL    4:       INCREASE    PRECIPITATION    CURING    MARCH    TFRCUCH   SEPTEMBER    3Y 
C***  THE    FOLLOWING    PERCENTAGES    -    75?    FOR    RAINFALLS    <    0. 10    IN, 

C**»  30?    FOR    RAINFALLS    <    0.50    IN,    10?    FCR    RAINFALLS    <    1.0    IN, 

C»**  NO    ACJUSTMENT    FCR    RAINFALLS    >    1.0    IN 

C 

0032  IFIMCDEL.EO.l J    GO    TO    40 

0003  IFIMCNTH.LT.3)    GO    TO    10 

0004  IF{VCNTH.GT.9)    GO    TC    13 

0005  A0J=1.75 

0016  IFIPRECIP.GT.0.10)    ADJ=1.30 

0C07  IF(t>KECIP.GT.0.50J     A0J=1.10 

000"  IF(WCCEL.E0.4)    GO    TO    20 

0C09  i'tPREClP.GT.1.0)     A0J=0.90 

0010  GO    TO    30 

0011  10  ACJ=1  .0 

0012  IF(MCDEL.EC.3J     AOJ=WPCNT 

0013  GO    TC    30 

0014  20  IF(PRECIP.GT.l.O)    AOJ*l.O 

0015  30  ?R£CIP=PR£CIP*ACJ 

0016  40  RETURN 

0017  END 
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FORTRAN    IV    G   LEVEL      21  SNGWRT                               DATE    =    79045                        21/27/40 

0001  SUBROUTINE    SNOWRT (PREC I ?, WATER .PACK , PET .TEMPAV  .SNCVAP) 
C 

C  *♦*    CALCULATION    OF    MOISTURE    ADDEO    TO    OISPCSAL    AREA    CUE    TO 

C  SNOWMELT    ON    THE    AREA    *** 

C 

0002  PEAL    M,MA,MR 

0003  M=0.0 

0C04  [F(PACK.GT.O.  1  )     SNQVAP=PET 

OC05  PACK=PAC<-SNCVAP 

0006  IF (SNCVAP. CT. 0.0)     PET=0.0 

0C07  IFJTEPPAV-32.)     10,10,20 

0008  10     IF(PPECIP)     70,70,30 

0009  20     IF(PACK)     90,90,40 

0010  30    PACK  =  PACK*PF<CCIP 

0011  WATEP=0.0 
0C12  on    TC    90 

C»*»      ma    !S    SNOWMELT    DUE    TC   ATMCSPHERIC    CONOITICNS 

0013  40    MA=0.C5*(TEMPAV-34.  ) 

0014  IF(MA.LT.O.O)     >*A=0.0 

0015  IF(PACK-MA)     60,60,50 

C«*«       MR    IS    SNOWMELT    SUE    TO    RAIN 

0016  50    MR=(PREClP*{TEMPAV-32.n/144 

0017  M=MR*MA 

0018  IF (PACK-*)    60,70,70 

0019  60    M=PACK 

0020  "ACK=0.0 

0021  GO    TC    80 

0022  70    PACK=PACK-M 

0023  SO    WATER  =  M«-?RECIP 

0024  90   RETURN 
0  02  5  END 
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FORTRAN    IV    G   LEVEL       21  OARCRT  DATE    =    79C45  2i/27/<«0 

0001  SUBROUT  INE    CARCRT IPERC . FCU.  SMUZ.FCL  ,  SMLZ  ,  FCCW  ,  SMGiiZ  ,DPERC,CCNCUZ  . C 
10>iOLZ  ,CGNDGVi,HUZ,HLZ,H3W,  PWU,  PWL,  PWG) 

0002  OnFNSICN    H(  3)  ,HK3)  ,SMO<  3)  ,WF{3)  .C0N013)  ,DE?(3},C12) 

0003  DATA    DEP.RCHGS,  EXCESS,  SCHGR/30.  48,  91.44,60.  96,  2*0.0/ 
C 

C***      DISTRIBUTION    OF    WATER    AOOEO    TC    EACH    PLOT 
C 

0034  EXCESS=0.0 

0005  IF (FFRC.LF.O.O)    GO    TO    10 

0006  SMAVLZ=FCU-SMUZ 

0007  IF(S.«AV'JZ.LT.G.G)     SMA\»UZ  =  G.O 
0003  EXCESS=PEF-.C-SMAVUZ 

0009  IFIEXCESS.LT. 0.0)    EXCESS=0.0 

OCIO  SMUZ»S'1UZ  +  PERC-EXCESS 

0011  SMAVLZ=FCL-S«LZ 

0012  IFISMAVLZ.LT. 0.03  SMAVLZ=0.0 

0013  EXTPA=EXCESS 

0014  EXC£SS=EXC£SS-SMAVLZ 

0015  IF(EXCESS.LT.O.C)    EXCESS=0.0 

0016  SMLZ=SMLZ+EXTRA-EXCESS 

0017  SMAVC=FCGW-SMGW 

0018  IFISHAVGW.LT. 0.0)  SMAVGW=C.O 

0019  EXTPA=EXCESS 

0020  EXCES3=EXCESS-SMAVGW 

0021  IFIEXCESS.LT. 0.0)    EXCESS=O.C 
0  022  SmG'*Z  =  SMGiZ  ♦EXTRA-EX  CESS 

0023  RCHGK=EXTRA-EXCESS 

0024  IF(PCHGR.LT.O.O)    RCHGR=0.0 
C 

C**«      MOISTURE    REDISTRIBUTION    USING    THE    ONE-0  IMENSICNAL    CARCY    EQUATION 

C***       FOR    UNSATURATED    FLOW 

C 

0025  10    LCCUNT=1 
0C26  DTI«£=0.1667 

0027  IF(PERC.LE.O.O)     0TIM£*1.0 

0023  IF (PE3C.LE.C.0)    LCOUNT-6 

0029  OPE?C=EXCESS 

0020  WF(1)=PWU/12. 

0031  WF(2)=PWL/36. 

00"*2  wF(3)=PwG/24. 

0033  SMO(  1)=SMUZ/12. 

0034  SMD(2)  =  S,wLZ/36. 
0C35  SMD(3)»S*GViZ/24. 

0036  DC    20    K=l,2 

0037  20    IF(SMG(K).GT.1.0)     SMC<<)=1. 
C 

0***      CALCULATE    SOIL    MOISTURE    TENSION    IN    CM 

0033  Hl(l)=SMO(  D-HUZ 

0039  HM2)=SM0(  21-HLZ 

0C40  HI(3)=SM0t3J-HGW 

0041  DO    30    K=l,3 

0042  HU)  =  £XP(Hi  (K.)  ) 

0043  IF(H(K) .GT.1500.)    H<K)=1500. 

0044  IFtH(K) .LT.O.O)    H<K)=0.0 

0045  XXX=SMO(K)-WFU) 

0046  IFIXXX.LT.O.OU     SMQ  1  K  )  =WF  (K  )  *0.01 
C 

C**»      CALCULATE    UNSATURATES    HYDRAULIC    CONDUCTIVITY    IN    CM    PER    DAY 
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FCRTRAN    IV    G   LEVEL       21  OARCRT                                 DATE    =    790*5                         21/27/40 

00'+ 7  IF(K.EQ.l)     CCNC(K)=CCNOUZ*(EXP<  72.339*SM0U)  J  ) /I  SMO<  .O-rfF  I  K)  ) 

0  046  IFU.EQ.2)     CCMDtKl  =  CCNDLZ«(EXP<75.595*SMC(iai  )  /  (S  MC  (  K  i-  W  F(  K.  )  J 

0  049  IFIK.EQ.3)     CCN01K)  =  CCN00H»1EXP170.588*SM0(KJ )) / IS *0 { X) -HF U )) 

0050  IFICOrOiK) .GT.IO.O)     CCN0(K)=10.0 

0051  30    IF(CCNDU)  .LT.l.OE-07)    CONO  (K  )  *  1  .OE-07 
C 

C*«*   CALCULATE  PCISTUAE  FLCWt  IN  INChES 

0052  40  00  50  K  =  l,2 

005  3  QU)=(CUN0lK>»CCN0lK*in/2»CTIME*(H(K*l.)-H<)<l*0EP<K))/DEPm 

0054  50  CU)=Q(K)/2.54 

0055  S*UZ=SfUZ-Q( I) 

0056  SMLZ=SWLZ*C(l) 

0057  SMLZ»S/»lZ-Q(21 

0058  S«GWZ=S^LZ+a:2) 
0C59  RCHGS=RCHGS*Q(2> 

0060  LCCUNT  =  LCCLNT«-1 

0061  IFILC0UNT.LE.6J    GO    TO    40 

0062  OPEKC=DPERC*RCHGS+RCHGR 

0063  RETURN 

0064  ENO 
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FORTRAN    IV    G  LEVEL       21  STCRAG                               CATE    *    7SC45                        21/27/40 

0001  SUBROUTINE    STCR AG(P 1 ,P2tP 3t PR  EC  IP , SNCW ,FRCZE ,MCNTH,GRCW  ,QCRM, RUNOF 
lF,h*lNPT,RUNCFT) 

C 

C 

C  «**  CALCULATION  OF  FEEDLOT  RUNCFF  »♦< 

C 

0002  IFIBYPASS.E0.3)    GO    TO    60 

C**«       CALCULATE    3    DAY    ANTECEDENT    MOISTURE 

0003  A1=PI+P2  +  P3 

0004  P1=P2 

0005  P2=P3 

0006  P3=P^EC!P 

0C07  IF  (SNOW. GT. 0.0. AND. FROZE. EQ.OJ  GO  TO  10 

0C03  IF(PPECIP.LE.O.C)  GO  TC  5C 

0009  IF  (FRCic.EO.l)  GO  TO  40 

0010  IFIAM.LE.0.5.AN0.PRECIP.LE.0.5)  GO  TO  50 

C***   CALCULATE  FEFDLOT  RUNOFF  USING  3  CAY  ANTECEDENT  MOISTURE  CONOITICNS 
C***   MCOIFICATICN  CF  THE  SCS  METH3C 

0011  10  AMl=AM*?RECIP 

0012  P*ESIP=PRECIP+SNCW 

0013  RC='97.0 

0014  IF(MGNTH.LT.4.CR.M0NTH.GT.10)  GC  TO  ZO 

0015  IFIAM.LT. 0.751  RC=9L.O 

0016  IF!,'.«l.GT.GPCk«. AND. PRtCIP.GT. GROW  )    PRES:P=GROW 
C017  GO    TC    30 

0018  20    IFUM.LT.0.50J     3C=91.0 

0019  IF    ( SNOW. GT. 0.0)    RC=97.0 

0020  IFIAMI.GT.DGRM.  AND. PRECIP.GT. DORM)     PRESIP=DCRM 

0021  30    CS=l0C0.0/SC-i0.0 

0022  RUMrFF=(PKESIP-C.2*CS)**2/( PR  ES  I  P«-0.3*CS ) 

0023  RU"l10FF  =  RUNCFF»PRECIP-PRES  l?*S.\CW 
0324  SNGU=0.0 

0025  IF(RUriCFF.GT.0.C6)     RUNCFF=RLNCFF-0. 06 

0026  IF(PRESIP-0.2»CS.LT.C.O)    GO    TO    50 

0027  GO    TG    60 

0029  40    SNCW=SNCW*PRECIP 

0029  50    RJNOFF=0.0 

003C  IFIBYPASS.60.2)    GO    TC    70 

0031  60  SUN0FT=*WINPT/13630»7.48) 

0032  •  70  RETURN 

0033  END 
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ABSTRACT 

Providing  a  dependable  supply  of  water  for  industrial,  domestic  and 
agricultural  use  in  the  United  States  has  increased  in  importance  during 
the  past  few  years.   One  method  which  is  gaining  in  acceptance  for 
increasing  the  water  supply  in  a  particular  locality  is  weather  modification. 
This  research  was  conducted  to  study  the  effects  of  a  weather  modification 
program  on  the  hydrological  environment  in  Kansas.   Ten  locations  were 
tested  using  two  rainfall  alteration  models  which  produced  variable 
percentage  changes  in  the  natural  daily  rainfall. 

In  addition  to  supplying  new  sources  of  water,  a  method  to  increase 
the  efficiency  of  available  water  supplies  was  investigated.   The  method 
of  optimizing  available  water  supplies  considered  in  this  study  was  the 
utilization  of  effluent  from  wastewater-treatment  plants  for  irrigation 
purposes.   The  scope  of  this  study  included  the  development  of  a  pro- 
cedure to  design  wastewater  control  facilities. 

A  multipurpose  computer  program  was  developed  to  evaluate  the  effects 
of  rainfall  augmentation  on  surface  water  hydrology  and  to  provide  an 
efficient  means  of  designing  a  disposal  control  facility.   The  computer 
program  is  capable  of  evaluating  the  hydrologic  response  of  various  crops 
and  soils  to  added  rainfall  and  of  considering  selective  land  disposal 
and  evaporative  techniques  to  optimize  a  disposal  system. 

The  rainfall  alteration  models  resulted  in  two  to  three  inch  increases 
in  annual  rainfall  throughout  the  State.   Most  of  this  added  rainfall  results 
in  increased  evapotranspiration,  with  smaller  increases  noted  in  surface 
runoff  and  deep  percolation. 


